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ABSTRACT 


In  this  report,  the  feasibility  of  obtaining  reliable  HF  communications  within 
and  outside  of  polar  regions,  especially  during  disturbed  conditions,  is  investi¬ 
gated.  To  establish  methods  by  which  propagation  modes  and  field  strength 
losses  can  be  found,  a  number  of  formulations  are  developed.  These  are 
evaluated  using  vertical-incidence  radio  measurements.  Also,  some  experi¬ 
mental  data  on  a  test  propagation  path  is  analyzed.  The  observed  values  are 
in  excellent  agreement  with  the  expected  ones. 

To  understand  the  behavior  of  the  ionosphere  during  both  quiet  and  disturbed 
conditions,  analysis  is  made  of  several  ionospheric  parameters.  Analyses  for 
a  large  number  of  disturbances  in  absorption,  sporadic  E,  and  F2-layer  critical 
frequencies  are  given.  Following  these  analyses,  data  on  HF  equipment  and 
systems  is  presented,  so  that  a  basis  for  selection  of  equipment  parameters  is 
available.  A  computer  program  is  then  developed  to  determine  the  optimum 
routing  of  traffic  after  the  traffic  load  and  available  channels  have  been  found. 
Some  examples  are  given. 

Finally,  a  completed  system  is  presented  for  which  the  maximum  possible  re¬ 
liability  is  found  in  an  example  situation.  A  traffic  load  is  assigned  to  a  hypo¬ 
thetical  network  during  the  height  of  a  polar  cap  absorption  event.  The  results 
show  that  communication  can  be  maintained  during  the  event,  but  at  a  reduced 
level. 
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I.  INTRODUCTION 


A.  STATEMENT  OF  THE  PROBLEM 

It  is  the  aim  of  this  study  to  determine  the  possible  wavs  by  which  communica¬ 
tion  to  and  within  polar  regions  may  be  maintained  thr  oughout  ionospherically 
disturbed  conditions.  At  present,  high  frequency,  1  to  30  Me,  radio  communi¬ 
cation  between  middle  and  high  latitude  transmitting  and  receiving  stations  is 
rendered  ineffective  during  ionospheric  storms  produced  by  intense  solar  dis¬ 
turbances.  The  partial  and,  at  many  times,  complete  blackout  of  commtmica- 
tion  links  is  caused  by  increased  ionization  produced  by  the  greatly  enhanced 
solar  X-ray  and  corpuscular  radiation  which  is  absorbed  in  the  D  region. 

In  most  instances  of  present-day  communication,  there  is  little  or  no  attempt 
to  avoid  absorbing  regions  or  to  utilize  anomalous  reflecting  regions  in  any 
systematic  or  knowledgable  way.  It  is  tacitly  assumed  that  all  communication 
links  are  simultaneously  affected  in  more  or  less  the  same  manner.  However, 
analysis  of  ionospheric  data  has  shown  that  the  spatial  locations  of  the  absorb¬ 
ing  and  reflecting  layers  vary  in  a  coherent,  but  complex  pattern  throughout 
the  lifetime  of  a  storm.  Furthermore,  the  analysis  suggests  that  the  patterns 
occur  in  a  predictable  manner. 

The  problem  then  is  to  determine  practical  ways  by  which  the  communication 
difficulties  associated  with  ionospheric  storms  can  be  overcome.  Another 
way  to  treat  the  problem  is  to  determine  the  optimum  way  in  which  the 
disturbed  ionosphere  should  be  used  for  a  given  system  of  transmitters  amd 
receivers.  Viewed  in  this  way  the  possibilities  for  improvement  of  the  system 
itself  may  be  found.  It  may  be  that  in  order  to  achieve  a  specified  reliability 
certain  alterations  of  the  system  may  have  to  be  made  before  the  optimum 
operating  procedures  can  give  tlte  desired  results. 


B.  SCOPE  AND  LIMITS  OF  THIS  REPORT 

This  study  will  include  factors  needed  to  solve  problems  related  io  HF  communi¬ 
cations,  and  especially,  those  problems  fou.nd  c'-j;rtng  ionospheric  storms.  At 
this  time,  the  studies  undertaken  are  in  progress  so  that  only  preliminary  results 
can  be  given.  However,  attempts  will  be  made  to  show  the  overall  system  in 
as  complete  a  form  as  possible.  It  is  not  intended  under  this  contract  to  go 
into  the  details  of  transmitting,  receiving,  or  data-processing  equipment 
characteristics,  but  only  the  general  system  parameters  are  considered. 

It  is  intended  that  in  the  next  report  prepared  under  this  contract,  revisions 
will  be  made  to  include  improved  and  additional  results.  The  authors  would 
appreciate  helpful  suggestions  and  criticisms  of  the  report. 
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IL  GENERAL  THEORY 


A.  RADIO  WAVE  PROPAGATION  MODES 

HF  communication  between  twodistant  stations  on  the  Earthis  possible  onlyif  the 
wave  is  reflected  from  the  ionosphere.  This  region  is  located  bolwesn  about 
90  and  500  km  above  the  Earth  and  is  a  weakly  ionized  gas.  It  is  divided  into 
layers  designated  D,  E,F1 ,  and  F2  in  order  of  increasing  height  as  shown  in 
fit  are  1. 

When  a  wave  passes  through  the  ionosphere  the  electrons  are  set  into  forced 
oscillation.  The  positive  ions,  because  of  their  large  inertia  are-practically 
uneffected.  The  coupling  of  the  wave  with  the  electron  motion,  affects  the 
phase  and  group  velocity  of  the  wave,  and  is  described  in  terms  of  the  refrac¬ 
tive  index.  If  collisional  damping  and  the  Earth’s  magnetic  field  are  neglected, 
the  index  of  refraction  is  given  by 


whore  the  phase  and  group  velocities  are  given  by 
c 

V-  •  —  «na  V,  ■  us 

“ft  ® 

where  c  s  velocity  of  light  in  free  space. 

Reflection  of  a  vertically  incident  wave  takes  place  when  the  refractive  index 
approaches  zero,  i.  e. ,  the  group  velocity  approaches  zero.  Thus,  if  the  max¬ 
imum  electron  density  of  an  ionospheric  layer  is  then  the  highest  fre¬ 

quency  which  can  be  reflected  is  given  directly  by  equation  (2). 


In  ionospheric  terminology  this  is  known  as  the  critical  frequency  of  the  layer. 

These  equations  are  useful  in  solving  propagation  problems  even  though  they 
are  not  Uieoretically  satisfying  since  two  important  properties  of  the  ionosnhe.  e 
have  been  omitted  in  the  derivation,  namely,  the  Earth's  magnetic  field,  and 
collisional  damping.  These  problems  are  treated  in  the  magneto-ionic  theory. 
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Simply  stated  the  magneto -ionic  theory  is  a  solution  of  the  wave  equation  in  the 
presence  oi-a  static  magnetic  field  which  introduces  anisotropy  into  the  index 
of  refraction  somewhat  analagous  to  a  birefringent  crystal.  Ai  a  result  the 
incident  wave  emerges  as  two  waves  elliptically  polarized,  one  clockwise  the 
other  counterclockwise.  As  in  physical  optics  they  are  tcrmcci  extraordinary 
waves  (fx)  and  ordinary  waves  (to)  respectively.  The  difference  in  their  crit¬ 
ical  frequency  is  given  approximately  by 


Other  results  of  the  magneto-ionic  theory  will  be  introduced  when  necessary. 
The  reader  is  referred  to  Budden  (1961)  or  Ratcliffe  (1959)  for  a  definitive 
treatment  of  the  problem. 

To  determine  the  behavior  of  a  wave  obliquely  incident  on  the  ionosphere, 
Snell's  Law  of  geometrical  optics  is  applied.  This  law  states  that  at  the  inter¬ 
face  of  two  homogeneous  layers  a  wave  passing  through  will  be  refracted  ac¬ 
cording  to  the  relation 


where  n  is  the  refractive  index,  <f>  is  the  angle  of  incidence  at  the  interface, 
and  the  subscripts  indicate  layers  1  and  2.  If  we  consider  the  ionosphere  to 
be  many  horizontally  stratified  layers  then  we  may  write 


Ml  •‘‘‘h  Mi-1  *'“^i 
M2  sin^j  ’  sin^j.j 

where  i  refers  to  the  ith  layer. 

By  successive  substitution 
Ml  *‘"^i 

>  TT  -  .  (5) 

Mn  i**2 

in  the  case  of  the  ionosphere,  o  is  taken  as  imiiiity  so  that  equation  (5)  refers 
to  a  continuously  varying  media. 

If  <l>i  is  the  angle  of  incidence  at  the  ionosphere,  then  (<1  =  1.  So 
sia6, 

a _ _  (6) 


where  corresponds  to  the  index  of  refraction  at  the  top  of  the  wave  path,  it 
will  be  denoted  as  n .  Also  at  the  level  where  the  wave  has  reached  maximum 
height  (travelling  horizontally),  =  90°  and 


»  sin^j  . 


In  other  words  the  refractive  index  at  the  point  of  reflection  is  given  b/  the  sine 
of  the  angle  of  incidence  at  the  ionosphere. 

Combining  this  equation  with  equation  (1) 

/ 

- —j  ■  sin^  “  (1  - 

so 

2 

—  •  cos  ^ 

(7) 

f  sec  <f> 


where  f  is  the  frequency  of  propagation  and  f,  is  the  plasma  frequency. 


Equation  (8)  is  called  tho  secant  law  for  oblique  transmission. 


The  angle  of  incidence  in  equation  (7)  for  propagation  over  a  curved  Earth,  can 
be  found  from  simple  geometric  considerations  and  is  expressed 


(8) 


where 


R  =  radius  of  the  Earth 


h,  =  virtual  height  of  reflection 
0  -  angle  of  transmission 


Tlie  angle  of  transmission  (sometimes  called  the  "take-off  angle")  written  in 
terms  of  the  height  of  reflection  and  the  great  circle  distance  d  from  trans¬ 
mitter  to  receiver  is 
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d 

2R 


(9) 


2R 


A  graph  showing  the  angle  of  incidence  as  a  function  of  distance  and  height  of 
reflection  is  given  in  figure  2.  Similarly,  the  trwsmission  angle  is  given  in 
f'^ure  3.  These  formulas  and  graphs  are  considered  standard  aind  can  be  found 
in  a  number  of  textbooks  dealing  with  HF  propagation. 

In  terms  of  the  maximum  range  o'  cbocrved  critical  frequencies,  and  reflection 
heights  of  the  various  layers  of  the  ionosphere,  table  I  has  been  prepared  in¬ 
dicating  possible  reflection  modes  over  a  curved  Earth,  Not  all  of  these  possible 
modes  apply  to  any  one  frequency,  location,  or  time,  but  they  serve  to  indicate 
the  ionospheric  propagation  modes  which  may  be  responsible  for  a  signal  at  a 
given  distance  from  the  transmitter.  For  example,  at  nighttime  there  is  no 
regular  E  layer  or  FI  layer.  Therefore  the  only  regular  mode  of  propagation 
is  via  the  F2  layer. 

For  oblique  incidence  the  maximum  usable  frequency  (mui)  is  found  in  the 
following  way.  With  the  aid  of  equations  (4)  and  (9)  the  secant  law  is  written  as 


(h,  +  2Riio^  —1 


r 


7  d  ' 

+  4R(Rth,)  lin^  — 
'  2R 


1/2 


where  =  plasma  frequency'. 

The  plasma  frequency,  which  is  a  function  of  the  height  is  written 


(10) 


f,  -  g(hp  (H) 

and  is  given  directly  on  an  ionogram.  (An  ionogram,  discussed  later  in  this 
section,  is  a  record  of  a  vertical-incidence  sweep- frequency  sounding  of  the 
ionosphere. ) 

The  problem  is  to  determine  the  maximum  value  off  in  equation  (11)  such  that 
equations  (10)  and  (11)  are  satisfied  simultaneously.  To  do  this,  plot  f,  versus 
b,  for  different  values  off  (f  s  f|  ,(2  ,f})  as  shown  in  figure  4.  Next  an  ionogram 
trace  is  superimposed  on  figure  4.  Obviously,  the  points  at  which  the  curves 
cross  satisfy  both  equations.  The  maximum  value  of  f  will  be  found  by  the  curve 
which  is  just  tangent  to  the  ionogram  trace.  In  this  particular  case  the  maximvun 
t  corresponds  to  f  =  f j  . 
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Figure  2  ANGLE  OF  INCIDENCE  VERSUS  REFLECTION  HEIGHT  AND 
GREAT  CIRCLE  PATH  LFNGTH 
63-3599 
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Figure  3  TRANSMISSION  ANGLE  VERSUS  REFLECTION  HEIGHT  AND 
GREAT  CIRCLE  PATH  LENGTH 
63-3598 
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Figure  4  GRAPHICAL  COMPUTATION  OF  MAXIMUM  USABLE  FREQUENCE 

63-3597 


-10- 


1.  Ionospheric  Data 

The  electrical  characteristics  of  the  ionosphere  are  determined  on  a  regular 
basis  (every  15  minutes)  by  over  100  stations  in  the  Northern  Hemisphere 
and  by  a  smaller  number  in  the  Southen  Hemisphere,  each  usi:ig  vertical 
incidence  sounders.  Measurements  are  normally  made  automatically  over 
a  frequency  range  between  about  1  and  25  Me,  the  interval  being  covered 
continuously  on  most  equipment.  The  power  output,  receiver  sensitivity, 
sweep-time,  and- antenna  efficiency  vary-considerably,  especially  between 
equipment  of  different  national  groups.  Results  obtained  are  rather  v'iform 
for  regular  layers,  but  there  are  some  exceptions  as  with  sporadic  E. 

Fundamentally,  the  ionosonde  method  consists  of  sending  a  short  pulse  of 
energy  toward  the  ionosphere  at  a  given  frequency,  and  recording  its 
return.  Elapsed  time  is  measured,  and  the  (virtUctl)  height  of  reflection 
is  determined  assuming  that  the  radio  wave  travels  with  the  speed  of  light. 
Since  this  is  not  'strictly  true  near  the  height  of  reflection,  the  indicated 
(virtual)  height  is  not  the  true  height.  In  certain  instances  the  signal  may 
not  be  returned.  There  are  two  causes;  (1)  the  absence  o.^a  reflecting 
layer  for  the  frequency  involved  and  (2)  the  strong  absorption  in  the  D  region. 
In  (1)  the  waves  escape  outward  through  the  ionosphere  (wave  frequency  is 
greater  than  the  critical  frequency)  and  is  a  measure  of  the  maximum 
electron  density  of  the  layer.  Absorption,  however,  is  a  reduction  in  signal 
strength  below  the  sensitivity  of  the  equipment.  The  ability  to-receive  an 
absorbed  signal  is  dependent'on  the  receiver  gain,  antenna  gain,  and  trans¬ 
mitted  power,  and  these  are  characteristics  of  the  individual  station 
involved. 

2.  lonograms 

An  idealized,  typical  ionogram  is  presented  in  figure  5.  The  ionospheric 
parameters  which  may  be  determined  easily  from  an  ionogram  are 

a.  E  region  maximum  electron  density  (foE) 

b.  FI  region  maximum  electron  density  (foFl) 

c.  F2  region  maximiun  electron  d‘;n8ity  (foF2). 

Ur-ing  certain  at  sumptions,  one  ccin  also  determine  from  fmin  the  integrated 
attenuation  occur.ng  in  the  ionosphere  up-to  the  reflection  height.  It  is 
possible  to  compute  the  true  height  of  reflection  from  the  virtual  height, 
but  this  is  not  normally  done,  and  requires  laborious  calculation  unless 
a  digital  machine  is  used.  At  present,  the  only  electron  density  versus 
true  height  profiles  available  on  a  routine  daily  basis  are  Puerto  Rico  in 
the  CRPL-F  series  reports.  Monthly  mean  .(quiet  day)  electron-density 
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Figure  5  IDEAUZED  lONOGRAM  (ORDINARY  COMPONENT) 
63-3596 


profiles  have  been  determined  for  a  number  of  Western  Hemisphere 
stations  and  ionospheric  true-height  maps  constructed  (Wright,  Wescott,  and 
Brown).  Refer  to  Wright,  Knecht,  andDavies(1957)forafurtAerdescription 
of  types;  to  Spread  F,  Penndorf  (1962);  and  to  the  sections  on  sporadic  E 
in  this  report.  Scaling  procedures  are  described  indetailin  Wright,  Knecht, 
and  Davies  (1957)  and  Piggott  and  Rawer  (1961)  and  are  most  useful  in 
interpreting  ionograms. 

B.  NONDEVIATIVE  ABSORPTION 


To  calculate  the  signal  loss  in  the  D  and  £  regions,  consider  the  expressic..  for 
nondeviative  absorption.  The  signal  loss,  expressed  in  decibels,  is  given  by 


/.ht 

I.  >  8.7  p  scc<f>  I  k  ds 

i 


where  p  =  number  of  penetrations  through  the  D  region  and  if>-  incidence  angle 
of  the  radio  wave  ahd 


k 

where 


f  Ni- 

^1,2  4. 


(13) 


h,  a  height  of  reflection 


e  =  charge  of  an  electron 

.N  a  electron-number  density 

V  -  elci  iron- collision  frequency 
m  ~  mass  of  an  electron 

c  a  velocity  of  light 

H  a  refractive  index  (in  the  nondeviative  region  p  =  1) 

~  gyrofrequency 

a  a  angle  between  the  direction  of  propagation  and  the  ma;^,netic  field 

f  a  wave  frequency  (where  the  ±  sign  :n  the  denominator  corresponds 
to  the  ordinary  and  extraordinary  component. s  respectively). 


Integration  of  equation  (12)  using  an  electron  density  profile  of  Houston  (1956) 
and  the  collision  frequency  distribution  of  Nicolet  (1953)  shows  that  v  in  the 
denominator  can  be  neglected.  With  this  approximation  equation  (121  may  be 
written  (for  the  ordinary  component) 
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4.35  p  stc<j> 
me 


Ni/ 


.  ds 


(14) 


where  =  fjj  cosa  , 


If  f>  r  the  time  being  only  the  case  of  one  pass  at  vertical  incidence  is  considered, 
it  will  be  possible  to  develop  eq\iation  (14)  into  a  more  useful  form.  With  these 
restrictions  equation  (14)  is 


t  .  f’  N.dh 

■"“'■•’"•'bra.® 


(15) 


When  the  reflection  height  is  above  the  absorbing  region,  the  product  Ni'  may  be 
treated  as  a  constant  and 

L  -  --L-  (16) 

«  .  fL)2 

where  A  is  a  parameter  which  may  be  found  from  absorption  measurements. 

(Of  course  for  an  independent  determination  of  A,  the  height  distribution  of  both 
N  and  V  must  be  known. ). 


When  the  reflection  height  is  within  the  absorbing  region,  Nv  cannot  be  treated 
as  a  constant  but  must  be  expressed  as  a  function  of  height.  From  the  analysis 
of  electron  density  and  collision-frequence  profiles  an  approximation  Ni'  - 
a  (h-  ,0))  written  where  a  is  a  constant,  and  h(|^  ,  qj  is  set  equal  to  65 

(km).  (Sac  example  in  figure  6.)  Thus  equation  (15)  may  be  written 


L 


..  (h-65)dh 


A' 

(TTfj/ 


(h,-65)2 


(17) 


The  reflection  height  in  equation  (17)  is  found  with  the  aid  of  the  height  distribu¬ 
tion  of  electron  density  and  the  relation  N  -=  1.  24  x  10^  f^.  Empirically,  it  is 
found  that  frequency  versus  heights  of  reflection  is  given  by  a  curve  of  the  type 
shown  in  figure  7. 
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DATE  OF  ABSORPTION 


Figure  6  ABSORPTION  COEFFICIENT  VERSUS  HEIGHT 
63-3595 
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Figure  7  FREQUENCY  VERSUS  HEIGHTS  OF  REFLECTION  IN  THE 

D  REGION 
63-3594 


The  upper  portion  of  this  curve  (figure  7)  can  be  thought  of  as  being  nearly 
parabolic  whose  equation  is  as  follows: 

(f,-foE)  -  -b'(h^-ho)2  (18) 


where 

(f  =  frequency  at  reflection 
foE  =  E  region  critical  frequency 
ho  =  height  at  the  critical  frequency 
hp,  =  the  height  indicated  by  tha  parabola. 
Solving  for 


(19) 


For  reflection's  below  h^  the  negative  vign  is  chosen, 
f,  -  0  and  h„  -  hj  and  substituting  in  equation  (19) 


hju  =  hfl  +  (hj  -  hg) 


Evaluating  b'  for 


(20) 


where  hj  is  a  constant  related  to  the  parabola  described  by  equation  (18).  This 
equation  represents  a  curve  which  cuts  the  h-axis  at  hj  (dotted  line  in  figure  7). 
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To  modify  equation  (18)  so  that  it  agrees  with  the  solid  curve,  some  hyperbolic 
function  must  be  added,  such  as 


where  the  first  factor  on  the  right  guarantees  that  as  the  frequency  approaches 
foE,  the  magnitude  of  h,  (height  contribution  of  the  hyperbolic  function)  approaches 
zero  which  is  what  is  required. 

The  equation  for  the  complete  curve  giving  the  height  of  reflection  H,  can  now  be 
written  as 


h,  -  h,u  + 

Using  the  condition,  at  f  ■  0,  h  »  65 
h,(f  •  0)  ■  hjy  +  •  65 

which  leads  to  c  •  hj  -  65,  and 
I  ton- I, (h,-65) 

\  foE  )  'tf*  I 

Adding  (20)  and  (22) 


hi  -  65 

h,  -  h„  +  (hi-h<,)E‘/2  -  e2 


where 

E 


Substituting  (23)  into  (17)  yields 


(22) 


(23) 


L 


A' 

(f  +  fL)‘ 


|bo-65  +  (hj-h„)E*/2_ 


ij-65 


2 


(24) 


Analysis  of  electron  density  profiles  shown  in  figure  8  (Kane  1962)  indicates 
that  hg  i  105  km.  Inserting  this  into  equation  (24),  multiplying  the  left-hand 
side  by  (40/40)^  and  noting  that  f  »  f,  at  vertical  incidence 


This  equation  is  similar  to  the  inverse-frequency  square  law  for  nondeviative 
absorption  except  that  it  is  modified  by  a  correction  factor  (term  in  the  braces). 
This  correction  factor  is  introduced  because  the  height  of  reflection  is  not  con¬ 
stant  with  respect  to  frequency.  As  a  resulti  for  frequencies  less  than  foE  it 
predicts  a  lower  absorption  than  a  simple  inverse  square  law. 

The  parameters  A  and  b  in  equation  (25)  are  determined  with  the  aid  of  experi¬ 
mental  data.  The  value  of  A  indicates  the  intensity  of  ionization  in  the  absorb¬ 
ing  region,  while  the  parameter  b  governs  the  shape  of  the  electron  density 
profile.  The  value  of  A  may  be  found  by  employing  frequencies  greater  than 
foE  and  measuring  the  absorption.  In  this  case  the  correction  term  is  unity  so 

L  .  ^ _  (26) 

if  f,  >  foE  . 


For  low  and  middle  latitudes  the  monthly  medians  of  absorption  can  be  expressed 


.  (27) 

(f,  + 

The  value  of  A"i8  found  from  measurements  of  at  Freiburg,  Germany  (lono- 
spharcn-Institut  Breisach,  1958)  for  a  sunspot  number  of  around  IbO.  Analysis 
of  this  data  shows  that  for  one  pass  (p  «  I)  A"-  500.  With  A-  500  cos®*5  y,  equa¬ 
tion  (25),  when  applied  to  monthly  median  data  can  be  written 

A"co8°*5x  ,„i/2 


Having  determined  A  the  parameter  b  is  found  by  solving  equation  (28)  f'>r  a 
range  of  b's  and  comparing  the  results  with  measured  values. 


Measurements  at  Freiburg  plotted  against  ( f,  -t-  f^)  are  shown  in  figures  9 
through  13.  Superimposed  on  these  curves  are  the  values  of  L  for  various  val¬ 
ues  of  b.  From  these  curves  the  best  value  of  b  is  taken. 
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Figure  9  OBSERVED  ABSORPTION  AND  COMPUTED  ABSORPTION  FOR 
VARIOUS  VALUES  OF  b,  JANUARY  1958 
63-3589 


Figure  10  OBSERVED  ABSORPTION  AND  COMPUTED  ABSORPTION  FOR 
VARIOUS  VALUES  OF  b,  MARCH  -  APRIL  1958 
63-3590 
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Figure  13  OBSERVED  ABSORPTION  AND  COMPUTED  ABSORPTION  FOR 
VARIOUS  VALUES  np  h.  AUGUST  -  SEPTEMBER  1958 
63-3593 
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Figure  14  MONTH  OF  THE  YEAR  (1958)  VERSUS  b 
63-3358 
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Seasonal  variations  of  b  are  apparent,  and  the  estimated  variations  of  b  for  the 
year  1958  are  plotted  in  figure  14.  Months  with  a  sunspot  number  greatly  dif¬ 
ferent  from  the  yearly  average  (180)  were  omitted.  It  was  found  that  b  is  re¬ 
lated  to  X  by  the  expression 

b  »  0.375  cos  X  •  (29) 


In  oblique  transmission,  the  distance  traversed  by  the  wave  in  the  D  region  of  the 
ionosphere  is  increased.  This  increase  in  the  path  length  through  the  abso?'Mng 
region  is  proportional  to  the  secant  of  the  angle  of  incidence.  Thus,  A  in  equa¬ 
tion  (25)  is  replaced  by  A  sec  i  .  The  reflection  height,  on  the  other  hand,  is 
lowered  since  it  is  computed  on  the  basis  of  the  secant  law. 


and  equation  (25)  becomes 

L  -  -  bE‘/2  „  /L^ 


if  f,  <  foE. 

When  f[  is  greater  than  foE  equation  (26)  is  modified  only  by  a  factor  see  ^  to 
take  into  account  the  longer  path  length.  Thus, 

L  -  ^  01) 

«r  ^ 

if  f,  >  (oE  . 

It  should  be  noted  that  the  value  of  fj^  may  be  different  in  oblique  propagation 
from  its  value  in  vertical  propagation. 

The  A  in  both  equations  (30)  and  (31)  is  determined  from  fmin  dtta  exactly  as 
described  before,  with  the  vertical  incidence  measurements  being  made  at  places 
where  the  wave  passes  through  the  D  region. 

These  equations  for  absorption  have  been  verified,  so  far,  only  by  comparison  of 
computed  versus  measured  values  of  absorption  at  vertical  incidence.;  ^he  re¬ 
sults  of  which  are  shown  in  figures  9  through  13.  It  is  also  desirable  to  compute 
electron  density  profiles  by  the  equation  giving  the  reflection  height  ind  compare 
the  results  with  profiles  measured  with  rocket-borne  instruments.  By  substitut¬ 
ing  the  expression  b  ■  (105  -  hj)/40  (equation  25)  into  equation  (23),  can  be 
written 


hj  *>  40 


1  -  bE 


1/2 


1  -  b 

Tpnr 


+  65. 


Replacing  the  frequencies  by  the  corresponding  electron  deneiiies  using  equation 
(2).  the  electron  density  as  a  fimction  of  height  is  obtained.  It  is  expressed  as 


hj  »  40 


1-bN**^^  -  -  —  N*^ 

1 


+  65 


(32) 


where 

N’ 


N^1/2  _  n1/2 

TiTi 


Equation  (32)  can  be  evaluated  by  finding  h,  for  various  values  of  N  which  are 
less  than  or  equal  to  Ng.  Such  computations  are  made  using  values  of  Ng  (or 
(oE)  and  b  appropriate  to  the  time  and  dates  of  the  measured  profiles  shown* in 
figure  IS  (HeilskiUi  et  al. ,  1960),  figure  16  (Pfister  and  Ulwich,  1961)  and 
figure  17  (Kane,  1962).  As  indicated  on  the  figures  the  computed  values  give 
excellent  agreement  with  the  measured  ones. 

C,  REFLECTION  COEFFICIENT  FOR  SPORADIC  E 


Regular  ionospheric  layers  (E,F)  are  assumed  to  be  totally  reflecting  below  their 
critical  frequency,  and  transparent  above  the  critical  frequency.  This  is  the 
mathematical  result  of  applying  a  full«wave  solution  of  Maxwell's  equations  to 
a  semi -infinite,  horizontally  stratified,  monotonically  increasing,  slowly  vary¬ 
ing  ionized  region  with  negligible  collisions  (Ratcliffe,  1959).  Below  the  criti¬ 
cal  frequency  all  the  power  is  returned  while  above  it  none  is  reflected.  The 
fact  is  ubservationally  confirmed  in  .nil  but  a  limited  number  of  special  situa¬ 
tions  in  which  the  above  referenced  solution  is  not  valid.  Specifically,  this 
involves  the  Z  trace  observed  on  high  latitude  ionograms  (Rydbeck,  1951),  For 
communications  the  Z  trace  is  not  of  practical  importance. 

In  the  case  of  sporadic  E,  the  space  gradient  of  electron  density  is  too  high  for 
normal  ray  theory  to  be  valid;  sporadic  E  is  partially  transparent  below  the 
plasma  frequency  and  partially  reflecting  above  it.  Tho  layer  is  found  in  the 
region  between  100  to  120  km.  Some  causes  which  have  been  suggested  are 
turbulence,  in  Whitehead  (1962),  wind  shear  in  Axford  (1962),  or  ionization  due 
to  corpuscular  radiation.  The  latter  is  likely  to  be  the  cause  of  most  high- 
latitude  Es,  while  the  former  two  are  thought  to  be  operative  at  lower  latituaes. 
The  thickness  of  the  ionized  region  may  be  of  the  order  of  several  wavelengths 
at  HF  (Seddoni  1961).  An  assumption  of  infinite  thickness  is  not  warranted,  and 
rigorous  solution  of  the  wave  equation  must  be  attempted.  The  reader  is  re¬ 
ferred  to  Rydbeck  (1943),  Epstein  (1930)^  Northover  (1962)  for  several  of  the 
solutions.  One  of  these,  due  to  Northover,  is  selectsd  as  offering  the  most 
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promise  for  high-latitude  phenomenon.  The  solution  is  one  'which  permits  the 
use  of  the  secant  law  for  propagation.  It  is  thus  applicable  to  vertical  incidence 
as  well  as  oblique  propagation. 

The  ionospheric  layer  is  taken  to  be  parabolic  in  shape  with  a  half  thickness  r 
and  maximum  electron  density  Solution  of  Maxwell's  equations  yields  a 
reflection  coefficient  which  is  graphed  in  figure  18.  These  cuvves  are  almost 
indistinguishable  from  the  family  of  curves  which  satisfies  the  equatiori 

R  -  (33) 

where 

R  =  ratio  of  incident  to  reflected  voltages. 

{  =  operating  frequency  in  Me. 


It  is  possible  using  the  observed  value  for  C  to  determine  the  layer  '.hick- 
ness  parameter,  and  then  from  figure  19  read  off  the  reflection  loss  at  the 
plasma  frequency.  Using  an  experimental  reflection  loss  versus  frequency 
graph  the  plasma  frequency  is  determined.  The  plasma  frequency  is  not.  how¬ 
ever,  necessary  for  propagation  problems.  If  '.lie  paremeters  C  and  A  are  deter¬ 
mined,  slant  path  propagation  may  be  predicted  using  the  secant-law  reduced 
frequency. 

^  f 

The  reflection  coefficient  is  determined  by  experimentally  using  ionosonde  film 
recordings  for  the  hourly  gain  run.  In  this  section,  D  -region  abosrption  is  neg¬ 
lected:  all  data  is  read  from  twilight  or  nighttime  film.  Any  change  in  the  foEs 
of  the  lowest  Es  trace  or  its  multiple  reflections  (if  present)  must  be  due  to  a 
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FiEure  18  REFLECTION  LOSS  VERSUS  A  AND 
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Figure  19  REFLECTION  LOSS,  C,  AND  hn 
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partial  reflection  loss.  While  the  absolute  value  of  ionosonde  gain  is  uncertain, 
the  gain  run  changes  are  accurately  known  and  are  20  db/step  for  the  Canadian 
LiG-17  and  Cossor  Sounders.  The  gain  run  steps  in  the  American  C-2,-3,  and 
-4  sounders  are  a  ftinction  of  the  receiver  medium-gain  setting,  which  is  ad¬ 
justed  automatically  daily  and  seasonally;  the  average  change  is  10  db/step. 
When  considering  the  reflection  coefficient,  the  computation  is  not  broken  up  on 
a  per  pass  basis.  This  reflection  loss  may  be  e;q>res'jed  in  cerms  of  the  iono¬ 
sonde  sensitivity,  ground  reflection  loss,  D-region  absorption,  and  additional 
space  loss  for  the  multiple  reflections. 


R  - 

where 

-  20  log  N  -  G(N- 1)  -  2  NL  +  <p)/li 

G, 

= 

ionosonde  gain  a  40  to  50  db 

N 

3 

number  of  reflections  from  layer 

G 

= 

ground  reflection  loss  a  5  db 

L 

S 

D-region  absorption  -  per  pass 

a 

ionosonde  gain  change  HI-MED-LO  (db) 

If  the  data  are  taken  at  night  or  twilight  and  fsE*  for  the  reflection  or  its  multi¬ 
ple  being  considered  is  greater  than  about  6  Me,  the  D-region  absorption  may 
be  neglected;  its  effect  is  most  noticeable  on  multiples  N  >  5. 

Wright  and  Gautier  (1960)  report  an  approximate  verification  of  the  secant  law 
for  Es  on  a  slant  path  in  middle  latitudes.  Preliminary  data  from  the 
Washington-Boston  path  indicate  that  both  the  secant  law  and  equation  (33)  are 
valid  (refer  to  section  11,  paragraph  G). 

D.  TRANSMITTERS,  RECEIVERS,  AND  ANTENNAS 

Numerical  caluclations  of  signel  strengths  are  based  on  transmitted  as  well  as 
ionospheric  parameters.  The  factors  influencing  transmitted  signal  strength  are 
transmitter  output  power,  and  antenna  gain,  from  which  the  losses  in  the  trans¬ 
mission  line  and  antenna  are  subtracted.  The  field-strength  voltage  E  is  related 
to  the  energy  density  P  by  the  intrinsic  impedance  of  free  space  Rg  by 

E^/R,  -  P 
or 

■  E  -  (377  P)‘/2  ^  (,/„)  (36) 

An  isotropic  radiator  in  free  space  transmits  power  equally  in  all  directions,  so 
that  the  power  is  spread  over  a  sphere  of  area 

A  “  4ird^  • 

d  s  distance  from  antenna  in  km  and  thus  the  power  density  P  (d)  per  unit  area 
is 
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Pt 

4nd^ 


wattsAm 


2 


(37) 


where  P,  is  the  transmitter  power  in  watts  radiated  by  an  isotropic  antenna. 
Combining  equations  (36)  and  (37) 


5.45  X  10 


,-3 


,1/2 


,  v/m 


(38) 


with  d  in  km. 


Further,  when  a  transmitting  antenna,  which  has  a  gain  g, ,  is  used,  the  voltage 
in  (iv/m  is  givenrby 


E 


5.45  X  13^  ^  . 


(38a) 


and 


20  log  E  »  74.8  +  10  log  Pj  g,  —  20  log  d  +  10  log  g^ 


(38b) 


(db  above  1  fiv/m ,  with  d  in  km) 

In  some  field* strength  calculations,  the  transmitted  power  is  reduced  by  1/3 
and  called  the  effective  radiated  power  (NBS  Circular  462).  This  yields  a  signal 
strength  of  3  x  lO^ftv/m  at  1  km  for  an  effective  radiated  power  of  Ikw. 

1.  Antennas 


The  characteristics  of  antennas  which  are  of  importance  in  communicationc 
are  the  gain  (relative  to  either  an  isotropic  antenna,  or  to  a  half-wave 
dipole),  the  radiation  pattern,  and  the  effective  area.  The  latter  is  useful 
for  receiving  antennas,  while  the  former  two  are  employed  to  describe  an¬ 
tennas  used  for  transmitting  and  receiving  systems. 

In  a  receiving  system,  it  is  important  to  know  the  efficiency  with  which  an 
antenna  can  abstract  energy  from  a  wave  incident  upon  it.  Particularly  for 
a  half-wave  dipole  antenna,  the  current  I  through  a  load  in  its  center  induced 
by  a  uniform  electric  intensity  EoCv/m)  parallel  to  the  antenna  is  (Schelku.' off 
and  Friis,  1952), 
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V  (z)  dz 


+V4 


(Z  +  73) 


where 

=  signal  in  db  above  ons  n»  at  tecriver 
d  -  distance  in  km 
fyj  =  frequency  in  Me 

G,  =  receiving  antenna  gain  relative  to  hedf-wavc  dipole  (do) 

G(  =  transmitting  antenna  gain  relative  to  half-wave  dipole,  (db). 
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2.  Receivers 


The  characteristics  of  receivers  vary  with  the  manufacturer,  but  an 
acceptable  HF  general-coverage  communications  receiver  has  a  sensitivity 
Of  0.  5  ftv  at  the  antenna  terminals  for  a  signal-to-noise  power  ratio  of 
10  decibels  in  a  6  kc  bandwidth.  Signal-to-noise  ratios  necessary  for  com¬ 
munication  with  a  specified  error  rate  are  different  for  various  information 
coding  methods.  Certain  standard  signal-to-noise  ratios  have  been  estab¬ 
lished  experimentally  by  Watt  et  al.  (1958)  and  are  given  in  table  II, 

The  noise  power  available  at  the  antenna  from  atmospherics,  cosmic, and 
man-made  noise  is  combined  with  the  antenna  thermal  noise,  and  an  effec¬ 
tive  noise  figure  is  defined  sc.  that 

p.  =  (,  ktb  (41) 


where 


p„  =  external  noise  power 

=  effective  noise  figure  at  antenna  terminals 
k  =  Boltuman's  constant  =  1.38  x  10’  Joules/deg  K 
t  =  temperature  in  degrees  K 
b  =  Effective  bandwidth  in  cps. 

Losses  occurring  in  the  transmission  line  and  antenna  are  expressed  as 
equivalent  noise  figure  and  while  the  receiver  noise  is  f,  .  The  effec¬ 
tive  noise  figure  referred  to  the  receiver  terminals  is 


f  =  f,  -  I  .  f,  fc  f, 


if  the  temperature  is  the  same  in  all  networks,  Friis  (1944),  and  C.itchlow 
et  al,  (1955).  If  p„  is  the  minimum  signal  power  -equired  to  maintain 
satisfactory  service  when  the  signal-to-noise  power  ratio  is  t ,  vnen 

p,„  =  rfktb  (43) 

or  in  decibels 

P„  =  R+  F  +  B-  204  (44) 


TABLE  II 


REQUIRED  SIGNAL-TO-NOISE  RATIOS>!- 


System 

wpm 

Noise 

Dynamic 

Range 

Carrier 

C/N  1  kc 

^•ms  carrier  to  rms 
noise  in  a  1-kc  band 

Type 

db* ** 

db 

10% 

1% 

0.1% 

A  CW  good  operator 

15 

Tt 

21 

S 

-  1 

1 

(4) 

B  CW  good  operator 

12 

Att 

68 

S 

-11.  5 

-  2 

(4) 

C  CW  lair  operator 

12 

A 

68 

S 

-  3 

- 

-- 

D  CW  good  operator 

12 

A 

69 

S 

-10 

-- 

E  CW  good  operator 

16 

A 

69 

S 

-  8 

(0) 

-- 

F  CW  good  operator 

20 

A 

69 

S 

-  6 

(1) 

-- 

G  FSK  ±  17  cps 

40 

T 

21 

S 

-2.3 

0 

1.6 

H  FSK  ±  50  cps 

60 

T 

21 

S 

0 

2 

3 

I  FSK  ±  25  cps 

60 

A 

68 

S 

0 

10 

17 

J  FSK  ±  50  cps 

60 

A 

69 

S 

3 

11 

15 

K  FSK  ±  425  cps 

60 

A 

50 

S 

10 

18 

21 

L  FSK  ±  425  cps 

60 

A 

40 

ionospheric 

fading 

12 

23 

34 

M  FM-FSK  576  TTY 
channels 

60 /ch 

T 

21 

tropospheric 

fading 

43 

48 

53 

N  FM-  36  voice  channels 

100/ch 

T 

21 

tropospheric 

fad'ng 

49 

80 

*  As  measured  by  Watt  et  al.  (1938) 

**  0.  001  to  90  percent  in  a  1-kc  band 

Numbers  in  parentheses  are  extrapolated  or  based  upon 
extrapolated  data 
t  thermedl  noise 
tt  atmospheric  noise 
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vhere 


Pm  =  10  log 
R  =10  log  r 
F  =10  log  t 
B  =10  log  b  . 
and  10  log  kt  =  -204,  if  t  = 

In  HF  communications  systems,  certain  approximations  may  be  made.  The 
losses  associated  with  the  antenna  circuit  and  transmission  line  are  small 
in  the  HF  band,  and  the  limiting  factor  on  signal  detectability  is  external 
noise  of  man-made,  cosmic,  or  atmospheric  origin.  The  receiver  noise 
figure  (f  is  less  than  that  due  to  external  noise  at  all  frequencies  below 
50  Me. 

In  equation  (44)  we  may  then  replace  F  by 
F,  »  10  log 

with  no  loss  in  accuracy  and  equation  (44)  becomes 
P„  »  R  +  Fg  +  B  -  204. 

We  may  use  equation  (44a)  to  calculate  the  necessary  antenna  gains  and 
transmitter  power  necessary  for  a  given  signal-to-noise  ratio,  or  using 
a  particular  system,  determine  what  information  transmission  modes  will 
give  the  lowest  error  rate. 

That  is 


P,  n  L^j.  -  G,  -  Gj  +  Fj,  +  R  +  B  -  204 


(45) 


where  P^  is  the  radiated  power  in  watts,  Lj  =  space,  reflection,  and  absorp¬ 
tion  losses. 


E.  RADIO  NOISE 


I.  Atmospheric  Radio  Noise 


Probably  the  greatest  limitation  to  the  useful  sensitivity  is  the  presence  of 
atmospheric  radio  noise.  The  source  of  this  noise  is  from  lightning  dis¬ 
charges,  but  it  can  be  propagated  thousands  cf  km  as  with  ordinary  trans¬ 
mitted  signals.  The  level  of  atmospheric  radio  noise  is  usually  deter¬ 
mined  by  the  strength  of  the  several  world  thunderstorm  centers.  How¬ 
ever,  when  electrical  storms  are  in  the  immediate  vicinity  of  a  receiver, 
the  local  source  will  predominate  and  will  sometimes  greatly  exceed  the 
normal  level.  In  middle  latitudes  the  occurrence  of  thunderstorms  is 
limited  to  about  1  day  in  10,  and  in  high  latitudes,  about  1  day  in  a  100. 
Furthermore,  the  duration  of  these  thunderstorms  is  only  about  2  hours. 
So,  most  of  the  time,  the  dominating  sources  are  the  equatorial  ones, 

e.g. ,  Southeast  Asia,  Africa  and  South  America.  For  this  reason  the 
monthly  medians  of  measured  noise  levels  indicate  the  effects  of  the 
major  thunderstorm  centers,  rather  than  local  ones. 

An  attempt  to  compute  the  noise  expected  at  middle  or  high  latitude  sta¬ 
tions  requires  a  consideration  of  many  factors.  The  important  ones  are 
listed  in  the  following: 

a.  Radiated  power  from  lightning  discharges 

b.  Rate  of  lightning  occurrence 

c.  Size  of  thunderstorm  area 

d.  Distance  from  source  to  receiver 

e.  Loss  of  power  from  absorption 

f.  Loss  of  power  from  ground  scatter 

g.  Frequencies  which  can  be  propagated  by  ionospheiiC  reflection. 

In  the  calculation  of  the  power  output  P  of  the  thunderstorm  center,  several 
approximations  are  made.  One  is  that  the  average  number  of  lightning  strokes 
per  1-hour  thunderstorm  is  taken  to  be  200  (Handbook  of  Geophysics,  1960). 
The  time  duration  of  a  single  stroke  is  about  30  it  sec  (Pierce,  I960),  so  the 
fraction  of  time  St/T  of  noise  output  per  thunderstorm  is  St/T  a  200  x  30  x 
10-6/3600  =  1.67  X  10-6. 

The  average  power  output  in  watts  of  a  single  lightning  stroke  for  a  1  -kc 
bandwidth  is  given  by  Pg  =  1.6  where  f  is  in  Me  (Handbook  of  Geophysics 


-37- 


1960),  So  the  power  output  per  thunderstorm  is  PgSt  or  1.67  x  10“^  x  1.6  x 
=  2,67  X  10*6  f-4.  The  total  power  output  per  thunderstorm  region  is  P  -nPgSt/T 
where  n  is  the  number  of  thunderstorms  simultaneously  occurring  in  the  re¬ 
gion.  Thus  P  =  2, 67  X  10*6  nf“4.  The  value  of  n  depends  upon  the  region 
considered,  the  season,  and  local  time  of  day. 

The  ionospheric  propagation  modes  are  now  considered  in  order  to  ascertain 
the  contribution  of  each  thunderstorm  area  to  the  noise  at  a  particular -loca¬ 
tion.  As  an  initial  approximation,  the  mode  with  the  smallest  number  of 
ionospheric  and  ground  reflections  will  be  least  attenuated  and  thus  be  domi¬ 
nant.  This  corresponds  to  the  lowest  angle  of  radiation.  At  distances  -greater 
than  3000  km,  the  approximation  is  made  that  the  ray  path  in  km  .  the  approxi¬ 
mation  is  made  that  the  ray  path  in  km  is  equal  to  the  distance  along  the 
Earth's  surface.  Further,  the  low-angle  radiation  is  focused  by  reflection 
between  a  curved  ionosphere  and  Earth  and  results  in  a  gain.  Gf  .  of  approxi¬ 
mately  8  db.  Applying  these  corrections  to  the  expression  for  the  field  strength 
previously  developed,  we  find 

20  log  E  -  74.8  4  G,  -  20  log  d'  10  log  P  /ir/m  (d'  >  3000  km) 


where 

G{  s  focusing  gain  s  8  db. 

d'  3  distance  on  surface  of  the  earth  in  km 

P  3  total  thunderstorm  power  output  in  watts /kc. 

Other  important  losses  are  absorption  in  the  D  region  and  ground  reflection 
loss.  The  signal  voltage  after  N  hops  have  been  reflected  at  the  ground  N-1 
times.  Thus 


E 


E 


where  g  is  the  ground  reflection  coefficient 


(47) 


20  log  E  a  20  (N  —  1)  log  g  +  20  log  E  . 


{47a) 


A  typical  value  of  g  is  0.63. 
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In  the  calculation  of  L  ,  the  absorption  is  found.for  each  hopi  so  that 
N 


(48) 


1 


where  N  is  the  number  of  hopst  and  using  the  results  of  the  previous  section. 


f  A  sec 


(49) 


where  rA  is  dependent  upon  the  state  of  the  D  region,  and  refers  to  the 
angle  of  incidence  to  the  D  region.  For  long  distance  propagation  by  the  E 
layer,  sec  <Iiq  =  5,  and  by  the  F2  layer,  sec  a  3,rA  may  be  found  with  the 
aid  of  fmin  data  just  as  with  normal  transmissions. 

Combining  these  loss  and  production  •equations 

20  log  £  -  105.4  +  Gf  -  20  log  d'  +  10  log  P 


rA  sec0p 

+  20(N-l)logg - 20  log!  .  (50) 

(f  +  fL)^ 


The  resultant  is  the  rms  noise  voltage  in  r.  1-kc  bandwidth  at  the  antenna 
terminals  of  a  half-wave  dipole.  Computation  of  the  atmospheric  noise  figure 
at  the  antenna  can  be  accomplished  using  equations  (41)  and  (50),  Expected 
values  of  median  radio  noise  levels  may  be  found  in  C.C.I.R.  Report  No.  65, 
(1957). 

2.  Cosmic  Noise 

Cosmic  noise  is  made  up  from  essentially  randomly  polarized  waves  of 
galactic  origin.  The  sources  are  not  uniform,  but  are  located  in  several 
major  centers  in  the  celestial  sphere.  The  noise  level  away  from  these 
centers  is  somewhat  reduced.  Cosmic  noise  reaches  the  Earth's  surface 
at  frequencies  greater  than  the  F2 -layer  critical  frequency  multiplied  by 
the  secant  of  the  incidence  angle.  This  noise  may  arrive  at  a  receiver' 
either  directly  from  outer  space  or  indirectly  by  first  arriving  at  the 
Earth's  surface  and  then  propagating  by  the  usual  ionospheric  modes. 
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Above  about  25  Me  the  major  contribution  to  radio  noise  is  caused  by  cosmic 
sources,  while  below  this  transition  frequency,  man-made  and  atmospheric 
noise  predominates.  However,  in  polar  regions  where'both  the  foF2  and 
atmospheric  noise  may  be  low,  cosmic  noise  may  be  the  limiting  factor  at 
frequencies  as  low  as  5  Me. 

3.  Man-made  Noise 

When  atmospheric  noise  is  low,  radio  noise  generated  by  electrical  appara¬ 
tus  such  as  ignition  systems  or  electric  drills  may  be  the  limiting  factor- in 
the  useful  sensitivity  of  a  receiver.  Radar  noise  may  originate  at  any 
point  where  there  is  a  sudden  change  in  voltage.  This  condition  is  most 
frequent  at  electrical  switches  oi  any  type.  Radio  noise  may  also  originate 
where  an  electrical  discharge  occurs  -  as  a  result  of  ionisation  of  air 
surrounding  a  conductor  around  which  there  is  an  excessive  voltage  gradient 
(Merriman,  1954).  A  list  of  important  radio  noise  sources  is  provided 
below: 

Power  and  distribution  lines 

High-voltage  equipment 

Electric  railway  systems 

Wire  communication  systems 

Commutator  motors  and  generators 

Thermostats 

Fluorescent  lights 

Sign  flashers 

Oil  burning  furnaces 

Electrostatic  discharge 

Small  generating  plants 

Industrial,  scientific,  and  medical  apparatus 

Arc  welding  equipment 

Elevators 

Vehicles. 
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F.  CALCULATION  OF  RECEIVED  SIGNAL  STRENGTH 


In  ordinary  HF  communication  it  is  desirable  to  compute  the  received  signal 
strength  from  some  distant  transmitter.  Such  a  calculation  is  basic  to  this 
study  and -what  follows  will  serve  as  an  example  showing  how  the  calculation 
is  made  if  certain  parameters  are  first  specified.  Those  parameters  involve 
the  following  factors: 

1.  Equipment  characteristics 

2.  Distance  from  transmitter  -  space  attenuation 

3.  Possible  modes 

4.  Absorption 

5.  Other  losses 

6.  Radio  noise. 

The  basic  equipment  characteristics  required  are  the  effective  radiated  power 
along  with  the  antenna  gain,  and  the  receiver  sensitivity  and  its  antenna  gain. 
Included  in  the  receiver  sensitivity  is  the  allowable  slgnal>to«receiver  noise 
ratio  for  the  type  of  communication  service  under  consideration. 

Under  normal  conditions  the  receiver  noise  level  is  replaced  by  the  external 
noise  level.  For  example,  in  the  absence  of  external  noise  a  receiver  sen> 
sitivity  of  1/jv  means  that  a  received  signal  of  Ijiv  would  give  a  signal-to- 
noise  ratio  of  say  10.  However,  in  the  presence  of  atmospheric  noise  at  a 
level  of  one  fiv  the  signal -to>noise  level  would  only  be  one.  To  maintain  a 
signal -to-noise  level  of  10,  the  received  signal  strength  would  have  to  be 
increased  to  3 /tv  i.e.,  10  db.  So  it  can  be  seen  that  the  maximum  amount  of 
allowable  signal  loss  between  the  transmitter  and  the  receiver  is  the  differ¬ 
ence  between  the  signal  strength  just  after  transmission  and  the  minimum 
allowable  signal  strength  after  reception  .  Of  the  factors  determining  the 
signal-strength  loss  between  the  transmitter  and  receiver,  the  letst  varia.ble 
of  these  factors  is  space  attenuation.  Although  this  varies  slightly  with  the 
mode  of  propagation  it  is  essentially  dependent  solely  upon  the  distance  of 
propagation. 

In  determining  the  maximum  frequency  of  propagation,  it  is  necessary  ■.  < 
consider  the  layer  critical  frequencies  and  their  height.  Also,  this  is  neces¬ 
sary  to  determine  the  amount  of  absorption.  In  other  words,  before  -osses 
can  be  determined  it  is  necessary  to  find  the  modes  possible.  When  this  is 
done,  it  is  possible  to  compute  the  loss  from  absorption  by  making  use  of 
equations’ (29)  or  (30). 
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Other  losses  such  as  ground  loss  and  reflection  losses  associated  with  sporadic 
E  also  shotild  be  taken  into  account.  If  the  total  of  the  losses  for  a  particular 
frequency  is  less  than  the  maximum  allowable  signal  loss>  then  that  frequency 
may  be  used  successfully. 

G.  EXPERUylENTAL  DATA 


Field-strength  data  is  obtained  from  recordings  at  Wilmington.  Massachusetts 
of  WWV  transmission  at  2.5.  '5.  10.  15.  and  20  Me.  Records  of  these  meas¬ 
urements  for  a  period  of  a  day  and  a  half  are  shown  in  figure  20,  Inspection 
of  these  graphs  shows  several  interesting  features.  These  are  as  follows; 

1.  Daytime  absorption  on  2,5  and  5  Me. 

2.  Enhanced  signal  level  on  all  frequencies  around  0000  LT.  12  Jan  1963 
in,  connection  with  sporadic  E  layer. 

3.  Large  changes  in  the  signal  strength  of  10  Me  during  the  daytime 
over  short-time  periods. 

Explanation  of  these  features  is  readily  understood  with  the  aid  of  "(-plots," 
i.e.,  graphs  showing  the  basic  ionospheric  parameters  from  vertical  inci¬ 
dence  soundings  at  15-minute  intervals.  (  -plots  for  11  and  12  Jan  1963  taken 
at  Ft,  Monmouth,  New  Jersey  are  shown  in  figure  21.  The  parameters  and 
values  appropriate  to  the  features  listed  above  are; 

a,  fmm  at  1200  LT.  11  and  12  Jan  1963  is  2.3  Me  on  each  day. 

Using  the  results  of  the  section  on  absorption  the  absorption  formula 
as  applied  to  ionosondes  is 


L 


I 


r  A 

y 

(fmin  +  (l,) 


(51) 


With  Lj  =  18,  r  =  0.75  =1.4  and  fmin  =  2,3  the  value  re  .ulting 

for  A  is  328. 

For  oblique  transmission  at  5  Me  the  absorption  formula  is  written 


L 


r  A  p  sec 
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with  r  =  0.66,  p  -  2,  8«c  ^  s  3.2,  fx.  =  0.6,  f  a  5,  and  the  calculated 
value  of  A  s  328,  absorption  at  S  Me  is  Lj  »  44  db.  This  is  about  the 
same  as  the  difference  shown  on  figure  20  between  the  day  and  night 
level  of  signal  strength  at  5  Me.  (If  absorption  at  night  is  negligible, 
then  this  difference  is  approximately  the  absorption  at  noon). 

b.  When  foEa  is  greater  than  6.7  Me  (tec  s  3.0)  reception  on 
20  Me  is  possible.  As  shown  on  the  I  -plot  as  is  the  case  between 
0000  and  0200  LT,  12  Jan  1963  -  the  same  time  as  the  20-Mc  recep¬ 
tion  of  WWV. 

c.  When  foF2  is  greater  than  about  6.7  Me  (sec  ■  1.6)  reception 
of  WWV  on  10  Me  is  high,  otherwise  it  is  abeent  unless  sporadic  E 
is  present.  On  the  11  Jan  1963,  (oF2  is  sufficiently  high  between 
1100  and  1545  LT  except  for  a  brief  period  around  1330.  On  12  Jan 
1963  {oF2  is  sufficiently  high  only  around  1045  and  again  at  1200  LT. 
WWV  reception  during  the  daytime  as  shown  on  figure  20  corresponds 
closely  to  these  intervals. 
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m.  IONOSPHERIC  PATTERNS 


A.  UNDISTURBED  PATTERNS 


I,  E  layer 


The  E  layer  is  a  regular  ionospheric  layer  in  the  sense  that  the  critical 
frequency  foE  is  directly  related  to  the  solar  zenith  angle  x-  When  ioniza¬ 
tion  is  caused  by  photo-ionizationi  and  the  only  loss  mechanism  is  cauf  '.d 
by  recombination,  the  variation  of  electron  density  may  be  stated 


dN 

dt 


I  -  aN^ 


(52) 


where 

I  =  intensity  of  ionizing  radiation 
a  =  recombination  coefficient 
N  =  number  density  of  electrons  . 

Equilibrium  condiHons  are  assumed,  which  yields 


I  “  oN^  . 

(53) 

N=  ^ 

(53a) 

That  is,  the  electron  density  at  a  given  height,  under  equilibrium  condi¬ 
tions,  is  proportional  to  the  square  root  of  the  ionizing  radiation  intensity. 
Consider  for  a  moment,  the  photon  radiation  from  the  sun.  The  light  rays 
have  an  intensity  Ij,  per  unit  cross  section.  Projecting  into  a  norizontal 
plane,  the  radiation  intensity  becomes 

I  =  IqCosx 


and  from  (53a) 

N  ~  (cosx)^^^ 


(55) 
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Thus  the  ordinary  wave  critical  frequency  foE  isi  by  equations  (2)  and  (55) 
foE  -  A  cos*/'*  X 

where  A  is  a  function  of  the  sunspot  number  only>  and  is  tne  H  region 
critical  frequency  for  X  -  0*  (subsolar  point).  Equation  (54)  is  based  on 
an  assumption  that  there  Is  no  loss  iddntensity  of  ■oiai'  radi&Uoii  ••  it 
passes  through  the  ionosphere  abdve  the  point  described  in  the  Ersfloa. 

This  is  particularly  important  near  sunrise  and  sunset  and  is  included  in 
a  theory  due  to  Chapman  (1931)  and  yields  the  same  result  for  foE  when 
the  density  varies  in  accordance  with  the  barometric  equation  for  an  iso> 
thermal  atmosphere 

N  -  Noexp(-z/H)  (57) 

where 

N  =  number  of  molecules 
X  =  height  in  km 
H  =  scale  height. 

Sunlight  is  not  monochromatic  nor  is  the  scale  height  constant.  These 
problems  have  been  treated  by  Chapman  (1939)  and  Shimazaki  (1959)  but 
yield  more  complicated  results. 

An  empirical  expression  has  been  constructed  which  represents  the  varia* 
tion  of  foE.  Monthly  median  data  from  the  CRPL-F  series  was  selected 
to  represent  low.  medium,  and  high  latitude  stations  and  was  sorted  accoid- 
ing  to  solar  zenith  angle  and  smoothed  sun  spot  (SSSN).  A  good  fit  (within 
i;  3  percent)  is  obtained  with 

foE  «  Eg  +  Sco«x  foe  X  ^ 
where 

Eg  =  2.  65  Me 

S  -  (0.55  +  0.0066  SSSN) 

A  graph  showing  these  results  is  given  in  figure  22, 

2.  Fl  layer. 


The  Fl  layer  is  solely  a  daylight  phenomenon;  and  it  is  prominent  in 
equatorial  latitudes  in  all  seasons,  and  in  all  latitudes  during  the  local 
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Figure  22  foE  VERSUS  SOLAR  ZENITH  ANGLE  AND  SSSN 
63-3580 
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summer  months.  On  the  ionosonde,  the  frequencies  near  the  critical  fre- 
quencyt  sviffer  group  retardation  causing  a  cusp  in  the  ionogram.  Rawer 
(1956)  gives  a  simple  formula  for  foFl, 

foFl  -  K  (FI)  cos°*2  X'  (59) 

where 

K(F1)  -  (4.3  +  0.011  R) Me  . 

A  detailed  analysis  of  CRPL-F  noon  median  values  indicates  that  a  seasonali 
as  well  as  solar  zenith  angle  dependence  exists.  Although  this  complication 
has  no  immediate  explanationt  it  is  probably  related  to  long-term  changes 
in  the  F2  region.  The  formulation  to  be  used  is 

(oFl  -  4.25  -  0.007  X 

+  j^O.017  +  0.0015  cos  -0.001  cos  ^  [l  -  10“^  x]  SSSN  (60), 

where 

N  =  month  number  (e.  g.  •  January  »  1) 

SSSN  -  smoothed  sunspot  number 
3,  F2  layer. 

The  behavior  of  (oF2  when  undisturbed  is  more  complicated  than  with  (oF)  , 
Undisturbed  foF2  is  the  condition  present  when  the  F2-layer  critical  fre¬ 
quency  is  not  greatly  different  than  that  indicated  by  the  monthly  medians, 
provided  that  the  individual  values  of  the  median  show  a  relatively  small 
standard  deviation.  Since  ioF2  is  closely  correlated  to  geomagnetic  dis¬ 
turbances,  it  is  also  possible  to  define  an  undisturbed  F2-layer  condition 
as  that  which  is  found  when  there  is  only  little  or  no  geomagnetic  disturb¬ 
ance  present.  Undisturbed  foF2  can  be  predicted  in  a  way  similar  to  foFl 
or  foE.  The  value  of  foF2  depends  upon  local  Ume,  solar  zenith  angle, 
the  geographic  coordinates,  the  season,  and  sunspot  number.  The  local 
time  dependence  is  somewhat  different  than  with  fpE  or  foFl .  At  sunrise 
there  is  usually  a  sharp  increase  in  foF2 ,  but  at  sunset  the  change  is  a 
gradual  one.  The  geographic  coordinates  must  be  considered  for  two  rea¬ 
sons.  Whereas  the  geographic  dependence  of  foE  and  foFl  can  be  expressed 
as  a  simple  function  of  the  solar  zenith  angle,  the  geographic  dependence  of 
foF2  may  be  quite  different  from  that  indicated  by  a  simple  solar  -zenith- 
angle  relationship.  Also,  geographic  coordinates  must  be  taken  into 
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account  because  foF2  values  vary  partly  in  accordance  with  the  Earth's 
magnetic  field.  These  variations  are  sometimes  referred  to  as  a  longitude 
effect.  For  prediction  of  foF2  monthly  medians  a  method  has  been  worked 
out  by  Jones  and  Gallet  (1962).  Alsoi  use  of  monthly  median  maps  is  quite 
suitable  to  determine  normal  frequencies  of  operation, 

4.  -Sporadic  E. 


Sporadic  E  cannot  be  classified  as  a  regular  ionospheric  phenomenon.  That 
is,  one  cannot  demonstrate  a  1-to-l  correspondence  between  Es  and  some 
physical  cause  --  e.g.,  the  solar  visual>emi8Sion  spectrum.  This  is  pos¬ 
sible  with  the  D,  E  ,  and  F  layers,  In  spite  of  the  fact  that  there  are  areas 
where  sporadic  E  is  common  in  high  geomagnetic  latitudes,  it  is  unusual 
for  Es  to  be  present  more  than  half  of  the  time  at  a  given  hourly  measure¬ 
ment  (e.g,,  foEs  >  5  Me). 

The  sorting  criteria  (i.  e. ,  foEs  >  5  Me)  is  established  to  eliminate  the 
possibility  of  the  regular  E  layer  obscuring  a  sporadic  E  layer  embedded 
in  it.  Further,  it  provides  a  suitable  lower  bound  on  the  electron  density 
of  the -layer  and  permits  comparison  with  the  results  of  other  workers, 
(Smith,,  1957). 

To  illustrate  the  phenomenon,  maps  for  July  1958  have  been  prepared  for 
stations  north  of  50  degrees  geomagnetic  latitude,  and  are  included  as 
figure  23  (a  to  h).  Although  the  maps  are  presented  every  3  hours,  a  map 
was  constructed  for  each  hour.  The  data  are  not  normalized,  but  qualifying 
symbols  are  added  to  indicate  whether  blackout  or  equipment  trouble  re¬ 
duced  the  number  of  data  days. 

Various  ionosondes  placed  at  an  identical  geographic  location  can  be  ex¬ 
pected  to  report  different  values  of  (oEs.  To  ascertain  the.actual  effect 
of  differing  equipment  sensitivity,  an  analysis  of  1  day  of  the  film  on 
hand  from  the  ionosondc,  at  Churchill,  Canada  was  made.  As  suggested 
in  the  IGY-URSI  operating  schedule,  hourly  ionograms  were  made  using 
three  different  gain  settings.  Fortunately,  the  relative  gain  change  is 
accurately  known.  For  the  LG- 17  sounder  (in  use  at  Churchil  )  the  difference 
between  successive  gain  settings  is  20  db,  making  the  total  variation  in  gain 
40  db.  It  is  unlikely  that  the  nominal  gain  at  other  stations  differs  from 
that  at  Churchill  by  more  than  ±  10  db. 

The  raw  data  from  a  "typical"  daily  scaling  (figure  24)  indicates  q.i.clita- 
tively  that  the  ionosonde  gain  enters  into  the  reported  value  of  foEs ,  but 
it  does  not  change  the  number  of  occurrences  of  foEs  >  5  Me  by  n.ore  than 
±  10  percent  with  respect  to  the  Churchill  medium  gain  as  a  base.  While 
perhaps  altering  the  shape  of  the  pattern  contours  slightly,  it  does  not  in 
any  way  invalidate  the  general  results  discussed  above. 
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Figure  23  NUMBER  OF  OCCURENCES  OF  foEs  >  5  Me 
63-3579 
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Figure  23  (Cont'd) 
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Figure  23  (Cont'd) 
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Figure  23  (Concl'd) 
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During  auroral  disturbances,  stratified  Es  may  be  observed  at  heights 
ranging  from  115  to  220  km.  This  consists  of  one  or  more  layers  or 
sporadic  E  which  are  not  multiple  reflections  of  the  lowest  layer.  It  may 
be  shown  that  these  reflections  come  from  tilted  ionospheric  layers  at 
about  110  km  and  are  slant^path  propagation.  The  minimum  observable 
frequency,  fmin  of  the  several  sporadic  E  layers  arc  related  by 


fmin2 


h'\  ,  /hA 

~)  (fjmni+  in).  -  'h  ~) 
‘'2  /  \"2  / 


(61) 


where 

frainj  =  (min  of  lowest  layer' 

(min2  =  fmin  of  any  other  Es  layer 
hj'  =  virtual  heigh  of  lowest  Es 

t 

h2  =  virtual  height  of  any  other  Es  layer 

fjj  =  local  gyrofrequency  (magnetic  field  assumed  vertical) . 

If  the  lowest  layer  is  above  120  km,  it  is  probable  that  it  too  is  a  slant 
reflection. 


Analysis  of  sporadic  E  data  as  described  in  the  foregoing  shows  that  on 
the  average,  there  appear  to  be  several  recurrent  patterns.  The  primary 
one  is  the  propensity  of  sporadic  E  in  the  Canadian  auroral  zone  and  in 
Arctic  regions.  Secondly,  a  distinct  traveling  region  is  present  which 
encompasses  the  region  extending  ±  45  degrees  from  local  midnight.  The 
coincidence  of  the  above  described  zones  produces  a  nighttime  maximum  of 
considerable  importance  in  Canada.  A  third  sporadic  E  regio:.  is  found 
over  Scandinavia  during  the  hours  6600  to  1200  UT,  and  while  small  is 
apparently  significant  (Maehlum,  1961).  The  Canadian  maximum  is  the 
only  stationary  pattern  which  remains  distinctly  observable  tltroughout 
the  day.  A  more  precise  description  of  this  pattern  places  the  maximum 
approximately  15  degrees  east  of  the  line  joining  Baker  Laxe  and  Churchill. 
The  northern  extent  is  probably  not  more  than  5  degrees  north  of  Baker 
Lake;  it  certainly  does  not  reach  Resolute  Bay.  To  the  south  it  reaches 
a  position  somewhere  between  Churchill  and  Winnepeg,  and  it  appsers 
to  extend  further  to  the  east  than  the  west.  Seasonally,  the  solstice  months 
are  more  important  than  the  equinoctal;  however,  this  applies  with 
certainty  to  IGY  data  only. 
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The  midnight  traveling  auroral'zone  peak  is,  as  its  name  implies,  asso¬ 
ciated  with  the  commonly  defined  auroral  zone.  It  appears  mo;<t  intense  in 
the  sector  180°K  through  Canada  to  45®E.  Certainly  this  result  is  some¬ 
what  influenced  by  the  density  and  location  of  tlie  reporting  stations. 

Scandinavia  appears  to  be  another  preferred  region  for  Es.  though  consider¬ 
ably  less  extensive  than  the  Canadian  maximum,  and  furthermore  not  all 
stations  report  a  definite  daytime  peak.  The  temporal  maximum  lies  some¬ 
where  between  0600  and  1200  UT.  The  pattern  is  not  distinctive  during  the 
winter  months, 

A  summertime  and  equinoctal  phenomenon  known  as  "Thule  type"  Et  reaches 
a  maximum  at  about  17  UT  over  a  major  part  of  the  polar  cap  north  of  75® 
geographic  latitude,  Penndorf  and  Coroniti  (1958),  It  appears  to  be  a  quiet- 
tirne  phenomenon,  but  may  be  present  though  obscured  when  absorption 
events  (PCA)  are  in  progress. 

5.  Layer  Height 

To  compute  the  various  modes  of  propagation  available,  it  is  necessary  to 
know  the  virtual  height  of  reflection  from  the  E,  FI,  and  F2  layers.  Analysis 
of  vertical  incidence  data  in  middle  latitudes  shows  that  the  height  of  the  E- 
layer  base  is  very  close  to  100  km.  For  the  FI  layer  the  base  height  is 
about  210  km  on  the  average,  while  the  height  corresponding  to  the  maxi¬ 
mum  usable  frequency  for  hops  of  2000  or  3000  km  is  about  260  km.  An 
overall  average  height  may  be  taken  as  235  km.  In  the  case  of  the  F2  layer 
the  variations  in  height  are  greater  than  those  of  the  Fl  layer.  The  base 
heights  and  the  heights  corresponding  to  hops  of  3000  to  4000  km  are  given 
below: 


h'F2 

vertical  incidence 

h'F2 

3  to  4000  km 

Summer:  Day 

400 

460 

Night 

285 

460 

Winter:  Day 

400 

460 

Night 

285 

460 

For  an  overall  average,  a  value  of  400  km  may  be  used.  There  is  a  cer¬ 
tain  condition,  under  which,  the  height  of  the  F2  layer  may  be  much  greater 
than  any  of  the  values  listed  above.  This  occurs  when  the  difference  between 
foF2  and  {oFl,  denoted  by  M,  is  small.  The  increase  in  virtual  height  is 
caused  by  the  increase  in  group  retardation  while  the  wave  is  passing  through 
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the  FI  region  near  its  critical  frequency.  The  empirical  relationship  be¬ 
tween  the  F2  virtual  height  and  Af  is  approximately 


h'F2  = 


555 

Ti^ 


if  Af  <  3  Me  . 

B.  DISTURBED  PATTERNS 
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1.  F2  Layer 

There  are  a  number  of  types  of  F2  layer  disturbances  which  are  of  impor¬ 
tance  to  communications.  Probably  the  most  important  of  these  is  the 
world-wide  decrease  in  the  F2  layer  critical  frequency  associated  with  geo¬ 
magnetic  storms.  Details  of  such  disturbances  have  been  described  by 
Obayashi  (1954).  The  decrease  begins  at  about  the  same  time  the  main 
phase  of  the  geomagnetic  storm  begins.  The  specific  location  and  intensity 
of  depressed  values  depend  upon  the  universal  time,  the  time  elapsed  from 
the  beginning  of  the  main  phase,  and  the  main-phase  intensity.  Usually, 
the  recovery  of  foF2  to  its  norrhal  pattern  takes  a  longer  time  than  does  the 
the  recovery  of  the  geomagnetic  field. 

An  example  of  the  storm-time  behavior  of  the  F2  layer  is  provided  in  fig¬ 
ure  25  (a  to  f)i  see  Hill  (I960).  Isopleths  of  critical  frequencies  are  shown 
in  intervals  of  3  Me.  Departures  from  the  monthly  median  values  are 
shown  by  dashed  curves,  also  in  intervals  of  3  Me.  As  shown  in  figure  25, 
the  departures  from  foF2  monthly  medians  increase  in  value  and  geographi¬ 
cal  extent  even  after  the  maximum  of  geomagnetic  disturbance,  i.  e.  ,  around 
1200  UT,  13. Sept  1957.  The  most  depressed  foF2  values  are  found  some¬ 
time  around  0000  UT  on  14  Sept  1957. 

Another  type  of  foF2  disturbance  occurs  in  polar  regions  during  equinoctal 
and  winter  months.  The  outstanding  feature  of  this  type  of  disturbance  is 
an  anomalous  increase  in  the  value  of  foF2.  Such  increases  have  already 
been  reported  (Kamiyama,  1962).  In  the  present  work  anoma'ous  foF2  in¬ 
creases  are  studied  with  the  aid  of  .synoptic  mapr,  so  that  the  geographical 
extent,  the  development  and  movement,  the  relationship  to  the  normal 
pattern  may  all  be  found. 

Such  disturbances  are  shown  in  figures  26  through  29.  These  are  I'pical 
ones  selected  out  of  some  28  separate  cases.  In  these  disturbances  the 
high  daytime  values  of  foF2  extend  into  ^polar  regions  and  sometimes  well 
across  into  the  nighttime  side  of  the  polar  cap.  The  width  of  this  region  of 
high  values  is  about  2000  km  on  the  average.  In  just  about  every  case  the 
maximum  value  in  the  anomalous  region  is  within  1  Me  or  so  of  the  maxi¬ 
mum  value  found  in  middle  latitudes.  The  duration  of  these  disturbances 
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[.  OOOOZ,  14  Sept  1957 


Figure  25  (Concl’d) 


b.  0900  UT,  10  Sept 

Figure  26  foF2  ANOMALIES  IN  THE  POLAR  CAP 
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b.  0900  UT,  14  Sept 


Figure  27  foF2  ANOMALIES  IN  THE  POLAR  CAP 
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Figure  28  foF2  ANOMALIES  IN  THE  POLAR  CAP 
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b.  2300  UT,  29  Dec 

Figure  29  foF2  ANOMALIES  IN  THE  POLAR  CAP 
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d.  0100  UT,  30  Dec 
Figure  29  (Concl'd) 


varies  from  about  1  hour  to  a 'half  day,  but  for  the  latter,  the  most  intense 
period  is  about  2  or  3  hours.  The  frequency  of  occurrence  of  these  dis¬ 
turbances  may  be  rather  high,  but  whether  or  not  they  occur  every  day-- 
between  fall  and  spring  --  cannot  yet  be  said.  Their  occurrcj.ee  is  certainly 
not  restricted  to  times  of  geomagnetic  storms.  A  summary  of  significant 
information  relating  to  the  four  disturbances  shown  in  the  figures  is  given 
in  table  III.  In  each  of  these  cases,  except  the  one  of  14  Sept  1957  the  geo¬ 
magnetic  field  was  rather  quiet.  (It  may  be  noted  that  the  maps  for  this 
case  follow  the  series  of  storm  maps,  described  earlier,  by  8  hours.) 

TABLE  III 


foF2  ANOMALIES  IN  THE  POLAR  CAP 


j  Date  (1957) 

Time  (UT) 

Maximum  (oF2  (Me) 

mm 

0800 

9 

3 

0900 

11 

3 

14  Sept 

0800 

10 

5 

14  Sept 

0900 

10 

5 

12  Oct 

1800 

12 

3 

12  Oct 

1900 

12 

2 

29  Dec 

2200 

13 

3 

29  Dec 

2300 

8 

3 

30  Dec 

0000 

8 

3 

30  Dec 

0100 

7 

3 

2.  Absorption 


Three  general  geographical  patterns  exist  in  which  anomalous  absorption 
is  most  often-found.  These  are  low-  and  middle-latitude  daytime  absorp¬ 
tion  (associated  with  solar  flare  X-ray  emissions),  polar-cap  absorption 
(associated  with  solar-flare  cosmic  rays),  and  auroral  zone  absorptio:; 
(associated  with  precipitavion  of  charged  particles  from  the  outer  magneto¬ 
sphere). 

When  X-rays  or  energetic  particles  reach  altitudes  of  the  D  region,  they 
produce  a  large  number  of  electrons  and  ions.  It  is,  of  course,  the  ptesence 
of  these  electrons  together  with  a  high-collision  frequency  that  results  in 
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absorption  of  radio  waves.  Besides  the  production  rate  there  are  two  other 
processes  governing  the  density  of  D-region  electrons.  These  are  the  re¬ 
combination  rate  and  the  attachment  rate.  During  hours  of  sunlight  the 
attachment  process  does  not  play  a  role,  because  the  photo-detachment  rate 
is  also  very  high.  However,  at  night,  attachment  becomes  important. 

A  rather  complete  review  of  studies  on  sudden  enhancements  in  low-  and 
middle -latitude  daytime  absorption,  or  SID  is  given  by  Warwick  (1963). 

Here,  the  SID  is  described  only  in  a  general  way.  A  SID  begins  at  about 
the  same  time  a  flare  becomes  visible.  Apparently  the  solar  radiation 
causing  the  disturbance  occurs  only  while  the  flare  is  in  progress.  As  soon 
as  the  flare  subsides,  so  does  the  absorption.  Almost  always,  the  onset  of 
absorption  is  very  rapid,  while  the  decay  is  somewhat  longer.  Although 
not  every  flare  produces  a  SID,  the  most  intense  ones  are  usually  found 
with  large  flares,  i.  e.,  class  3  or  3+  .  However,  some  large  flares  do  not 
produce  very  much  absorption.  On  the  other  hand,  small  or  medium  size 
flares  can  produce  moderate  or  even  intense  bursts  of  X-rays,  which  lead 
to  corresponding  amounts  of  of  absorption.  The  duration  of  a  SID  varies  over 
a  range  of  time  from  a  few  minutes  to  2  or  3  hours.  An  example  of  a  mod¬ 
erate  SID,  as  represented  by  the  parameter  fmin  ,  may  be  found  on  the  chart 
for  31  Aug  1957,  0600  UT  shown  in  figure  30,  An  example  of  a  strong  one 
is  found  on  the  chart  for  18  Sept  1957,  1800  UT  shown  in  figure  31,  For 
communication  a  SID  may  cause  outage  on  links  which  are  in  daylight.  The 
lowest  frequencies  on  these  links  would  suffer  the  most  absorption.  The 
duration  of  such  outage  would  last  only  an  hour  or  two  at  most. 

Polar  cap  absorption  (PCA)  is  associated  with  solar  flares  but  in  a  different 
way  from  a  SID.  Following  the  occurrence  of  certain  flares,  which  can 
be  identified  by  the  flare  intensity  and  the  associated  solar  radio  noiiie, 
protons  in  the  1-mev  to  l-bi:v  energy  range  arrive  near  the  Earth  and  are 
guided  into  polar  regions  of  both  hemispheres  by  the  Earth's  magnetic  field 
(Obayashi  and  Hakura,  1960)  and  (Hakura  and  Coh,  1961).  The  sun-£arth 
travel  time  for  these  protons  is  about  2  to  6  hours,  but  in  some  instances 
the  travel  time  is  well  out  of  this  range  (Obayashi  and  Hakura,  1960).  The 
onset  of  polar-cap  absorption  falls  into  two  classes,  one  charseterized  by 
a  rapid  increase  in  absorption,  and  the  other,  a  slow  increase. 

The  geographical  distribution  of  polar-cap  absorption  for  typical  events  is 
shown  in  figures  32  and  33.  These  figures  have  been  taken  from  a  series 
of  synoptic  maps  of  eleven  periods  of  polar-cap  absorption  (Coroniti.  1963), 
The  solid  curves  on  the  maps  refer  to  values  of  fmia  in  units  of  Me.  The 
shaded  areas  indicate  regions  where  fmin  exceeds  foF2  .  This  condition 
means  that  there  are  no  signals  returned  to  the  ionosonde  in  sufficient 
strength  to  be  recorded.  The  dashed  curve  on  a  map  indicates  the  position 
of  sunrise  and  sunset  at  the  Earth's  surface.  Considering  the  patterns  of 
polar-cap  absorption  on  the  basis  of  where  the  energetic  protons  reach  the 
D  region,  it  would  be  expected  that  absorption  would  be  greatest  near  the 
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Figure  31  MAPS  FOR  SUDDEN  IONOSPHERIC  DISTURBANCE 
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Figure  32  (CONCL'IS) 


Figure  33  MAPS  FOR  POLAR  CAP  ABSORPTION  EVENTS 


Figure  34  MAPS  FOR  AURORAL  ZONE  ABSORPTION 


geomagnetic  pole  and  would  decrease  slowly  to  a  lower  geomagnetic  lati¬ 
tude  (about  65  degrees)  and  then  decrease  rapidly  further  away  from  the 
pole.  Such  a  pattern  is  actually  found,  but  with  one  important  exception, 
that  is,  on  the  darkened  portion  of  the  polar  cap.  Because  of  a  high  at¬ 
tachment  rate  the  free-electron  density  is  considerably  reduced.  Thus, 
a  solar-proton  event  of  a  given  strength  producer,  in  the  summeri  a  region 
of  absorption  covering  the  whole  polar  cap,  whereas  in  the  winter  the  region 
of  absorption  is  confined  mainly  to  the  relatively  small  illuminated  portion 
of  the  polar  cap.  Additional  information  on  the  extensive  work  performed 
in  the  last  few  years  related  to  PCA's  may  be  obtained  from  the  summary 
articles  by  Hultqvist  (1962)  and  Kahle  (1962). 

Toward  the  end  or  just  after  a  polar-cap  absorption  event  there  is  usually 
a  geomagnetic  storm  (Hill,  1961).  During  this  time  enhanced  absorption  is 
found  in  the  auroral  zones.  The  period  of  greatest  absorption  occurs  about 
the  same  time  as  the  maximum  of  the  main-phase  geomagnetic  disturbance. 
The  region  of  absorption  in  the  auroral  zone  usually  progresses  a  few  hun¬ 
dred  kilometers  in  a  direction  away  from  the  pole  (Obayashi,  1959).  In 
severe  geomagnetic  storms  this  displacement  may  be  as  much  as  2000  km 
(Hakura,  1961).  The  effects  of  attachment  are  not  so  great  as  with  polar- 
cap  absorption,  because  the  D-region  electrons  are  usually  produced  at  a 
higher  altitude  in  auroral-zone  events  than  in  PCA  events.  In  other  words, 
differences  in  ar.-iounts  of  day  or  nighttime  absorption  will  not  be  as  great 
as  with  polar  cap  events.  A  good  example  of  auroral  zone  absorption  as¬ 
sociated  with  a  geomagentic  storm  is  found, on  the  map  for  13  Sept  1957, 
0600  UT,  shown  in  figure  34. 

Auroral-zone  absorption  is  not  limited  to  periods  of  geomag.iCtic  storms 
of  a  classical  type,  (preceded  by  a  sudden  commencement  of  geomagnetic 
activity  and  a  polar-cap  absorption  event).  During,  periods  of  minor  geo¬ 
magnetic  activity,  which  occur  several  times  more  often  than  geomagnetic 
storms,  there  is  usually  a  portion  of  the  auroral  zone  where  there  is  in¬ 
creased  absorption. 

3.  Disturbed  Es 


Observed  values  of  sporadic  E  critical  frequency  during  geomagnetically 
disturbed  periods  (see  table  IV)  have  been  plotted  on  polar  stereographic  * 
maps.  The  maps  are  plotted  on  either  a  1  or  3  hour  UT  basis  and  isopleths 
of  foEs  drawn.  Three  latitudinal  zones  are  established  (U,  A,  L)  roughly 
8  degrees  wide  with  the  A  zone  centered  in  the  mean  observed  auroral  zone, 
and  the  U  and  L  zones  above  and  below  A,  respectively,  as  shown  in  figure 
35.  Longitudinally  the  zones  are  30  degrees  wide,  commencing  with  a 
zone  centered  on  the  Greenwich  meridian.  The  intensity  of  the  Esis  indi¬ 
cated  by  a  number  system: 
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Figure  35a  ZONE  DEFINITION  FOR  SPORADIC-E  ANALYSIS 

SPORADIC  E  CLASSIFICATION 

INDEX 

0 

1 

\% 

2 
3 

Figure  35b  CLASSIFICATION  OF  SPORADIC  E  ZONES 


RANGE  of  foEs  (MQ^) 


0*  foEt  •  4. 0 
4.0*foE$  *7.0 
4.0-foEs  *  10.0 
7.0*  foes*  10.0 
IO.O*foEs 
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Index 


Range  of  foEs  (Me) 


0  0.<  foEs  <  4.  0 

1  4.  0  <  foEs  <  7.0 

1-1/2*  4,0  <  foEs  <  10.0 

2  7.  0  <  foEs  <  10.0 

3  10.  0  <  foEs 


*The  value  1-1/2- occurs  in  prediction  only  and  results  from 
the  addition  of  the  spatial  and  auroral  patterns.  In  range  it 
it  encompasses  indices  I  and  2,  but  the  maximum  probability 
is  at  7  Me. 

From  the  maps  a  table  is  derived  which  can  conveniently  represent  a  series 
of  maps;  for  e^mple,  the  map  series  of  4  and  5  September  is  represented 
in  table  V. 

Each  map  is  reduced  to  a  single  line  in  this  table.  A  combination  of  these 
indices  listed  above  represents  the  intensity  of  Es  in  the  regions  of  U,  A, 

L,  within  a  given  longitudinal  zone,  e.  g. ,  U,  A,  L  equals  1,  0,  3. 

A  study  of  disturbed  periods  (listed  in  table  IV)  indicates  the  presence  of 
an  auroral-zone  traveling  of  maximum  Es  with  its  center  around  local  mid¬ 
night.  Its  latitudinal  and  longitudinal  breadth  are  similar  to  the  monthly 
patterns  except  that  during  the  main  phase  of  "great"  geomagnetic  storms, 
the  pattern  intensifies,  mov-;s  to  the  south,  and  spreads  over  a  greater 
portion  of  the  night  hemisphere.  This  is  also  noted  in  studies  of  the  polar 
aurora  (Akasofu  and  Chapman  1962;  Akasofu  1962).  In  terms  of  a  spatial 
pattern,  the  ionospheric  stations  at  Baker  Lake  and  Churchill  report  the 
most  persistent  sporadic  E  occurrence.  Sporadic  E  during  disturbed  pei'i'ods 
is  related  to  the  average  patterns  previously  described  (section  III.  A.  4. ) 
since  the  average  patterns  include  both  quiet  and  disturbed  intervals. 

Using  these  notions  a  prediction  table  is  constructed  for  sporadic  E  during 
disturbed  conditions  (table  VI).  This  prediction  table  consists  of  a  summa¬ 
tion  of  the  spatial  pattern  and  an  auroral-zone  traveling  maximum  Es  cen¬ 
tered  about  local  midnight.  (A  predicted  value  of  0.  5  in  the  total  pattern 
is  counted  as  0. )  This  sum  in  table  VI  iu  the  basic  predictionin  universal 
time.  There  are  also  storm-time  variations  which  are  dependent  on  the 
phase  of  the  storm.  Changes  in  numerical-value  assignments  for  both  the 
spatial  and  auroral  patterns  are  made  to  conform  with  the  storm-time 
changes  (table  VII). 
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TABLE  IV 


STUDY  OF  DISTURBED  PERIODS  (SPORADIC  E) 


Begin 

End 

Type  of  Map 

Date 

Time  (UT) 

Date 

Time  (UT) 

(Interval) 

19f ; 

3  Aug 

12 

4  Aug 

15 

3 

13  Aug 

06 

13  Aug 

10 

1 

3  Sept 

00 

3  Sept 

21 

3 

4  Sept< 

M 

5  Sept 

12 

1 

6  Sept 

00 

6  Sept 

21 

3 

13  Sept 

00 

13  Sept 

16 

1 

23  Sept 

00 

24  Sept 

06 

3 

6  Nov 

18 

7  Nov 

08 

1 

1958 

18  Jan 

00 

18  Jan 

21 

3 

11  Feb 

00 

11  Feb 

21 

3 

19  Feb 

00 

19  Feb 

21 

3 

11  Mar 

12 

12  Mar 

09 

3 

15  Mar 

00 

15  Mar 

21 

3 

22  Mar 

00 

22  Mar 

2i 

3 

17  Apr 

00 

17  Apr 

21 

3 

Time  UT  0  30  E  60  £  90  C  120  E  I  ISO  E  I  18)  I  ISO  W 
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TABLE  VU 

STORM  TIME  CHANGES  IN  THE  PREDICTION  SCHEME  FOR  Es 


Storm  Type 

Value  in 
table  VI* 

Minor 

Moderate 

Great 

Initial 

Phase 

Main 

Phase 

Initial 

Phase 

Main 

Phase 

Spatial 

0,  01  0.^0.  0 

... 

-- 

... 

... 

0.  0.  Si  0.  5 

— 

...  ' 

0.  5  0,  Ek  0.  0 

... 

... 

0.  5k  1.  a  0.5 

... 

0.  5  1.  ft  0.0 

... 

0.  ft  1.  ft  1.0 

0.  SkO.  SkO.O 

... 

... 

0.50.50.0  < 

... 

(0.  &  0.  5  0.  5) 

... 

... 

(0.  5'0.  5  0.  0) 

0.  ft  1.  ft  0.0 

O.ft  l.ft  1.0 

i.ai.ao.5 

... 

... 

2.  ft  1.  ft  0.5 

2.  ft  T.  ft  0.5 

1.  ft  1.  ft  2.0 

Auroral 

0.  a  1.  ft  0.5 

... 

... 

... 

... 

O.ftO  ft  1.0 

0.  5k  1.  ft  0.5 

... 

... 

... 

... 

O.ft  l.ft  1.0 

<<The  Values  in  this  column  change  to  the  values  indicated  in  the  other 
columns,  depending  on  stor^^  type. 

’)”!<Dashes  indicate  no  change. 
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As  an  example,  several  hours  selected  from  a  moderate  sudden-commence¬ 
ment  geomagnetic  storm  on  6  and  7  November  1957  are  presented  in  figures 
36  through  38  in  map  form  and  compared  with  the  predicted  values.  The 
storm  began  at  1821  UT  on  6  November  1957  and.had  a  sudden  commencement 
of  the  type  designated,  -SC-.  The  epoch  of  maximum  main-phase  de¬ 
pression  of  the  geomagnetic  field  occurred  at  0100  UT  on  7  November  with 
a  Dst  component  of  about  100  gammas. 

To  obtain  a  preliminary  indication  of  the  utility  of  the  prediction  scheme, 
the  total  is  found  of  the  number  of  times  the  prediction  does  not  coincide 
with  the  observed  value.  No  distinction  is  made  between  the  indices  1 
through  3.  There  is  a  total  of  36  prediction, zones. 


Time 

6/7  November  1957 

Number  of  Errors 

Percent  of  Error 

2100 

5 

14.0 

0000 

3 

8.3 

0300 

2 

5.6 

C.  RELATED  SOLAR-GEOPHYSICAL  DATA 
1.  Solar  Data 

The  intimate  connection  between  the  solar  disturbance  phenomenon  and  ion¬ 
ospheric  conditions  has  been  introduced  in  previous  sections.  Solar  sun¬ 
spots,  flares,  and  associated  radio  and  visual  emissions  create  or  influence 
disturbances  of  the  geomagnetic  field  and  ionosphere.  Parameters  which 
characterize  the  ionosphere  electron  density  have  been  previously  intro¬ 
duced.  Some  measures  of  solar  and  geomagnetic  disturbances  are  now  dis¬ 
cussed  which  are  used  in  correlation  studies  to  establish  the  mechanism  of 
the  interaction  between  the  sun  and  Earth  which  produces  an  ionospheric 
storm. 

Logically,  the  assumed  source  of  the  disturbance,  the  sun  would  first  be 
described  but,  in  fact,  it  is  beyond  the  scope  of  this  report  to  discuss  any 
of  the  theories  of  the  solar  reaction  mechanism.  Visual  and  radio  observa¬ 
tions  of  the  solar  disc  provide  a  measure  of  the  solar  disturbance.  The 
relative  sunspot  number  is  characteristic  of  the  total  activity  of  the  sun 
and  is  defined  as 

Rj  =  K(I0g  +  s)  (63) 
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Figure  36  OBSERVED  AND  PREDICTED  PATTERNS  OF  SPORADIC  E 
DURING  A  GEOMAGNETIC  STORM,  2100  UT 
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Figure  37  OBSERVED  AND  PREDICTED  PATTERNS  OF  SPORADIC 
DURING  A  GEOMAGNETIC  STORM,  0000  UT 


Figure  38  OBSERVED  AND  PREDICTED  PATTERNS  OF  SPORADIC 
DURING  A  GEOMAGNETIC  STORM,  0300  UT 
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where  g  =  number  of  sunspot  groups 

s  =  total  number  of  distinct  spots 

K  =  scale  factor  depending  on  observer. 

Two  different  relative  numbers  are  computed,  and  the  American  and 
Zurich  relative  sunspot  numbers.  The  magnitude  of  each  is  different,  but 
the  trend  is  the  same.  It  indicates  that  the  activity  of  the  sun  is  periodic 
in  about  11  years  with  a  probability  that  alternate  cycles  will  have  the 
close sf'amplitude  relationship.  A  12-month  running-mean  Zurich  smoothed 
sunspot  number  is  often  used  for  computational  purposes  and  termed  the 
smoothed-observed  sunspot  number  R. 

At  times,  small  regions  of  the  chromosphere,  almost  invariably  associated 
with  sunspot  groups,  suddenly  intensify  their  emissions  of  certain  spectral 
lines,  primarily  the  calcium  H  and  K  lines  and  the  hydrogen  Hg  line.  These 
flares  have  a  lifetime  of  a  few  minutes  to  hours,  with  the  average  being 
about  a  half  hour.  The  importance  of  flares  is  designated:.  1-,  1,  2,  3,  3+ 
primarily  upon  their  visual  characteristics.  There  is  a  good  correlation 
between  solar  flares  and  sudden  ionospheric  distrubances  (SID);  so  that 
these  ionospheric  absorption  events  can  be  used  to  identify  a  solar  flare. 
Normally,  this  radiation  effect  produces  abnormally  increased  HF  absorp¬ 
tion  and  is  the  cause  of  the  sudden  short-wave  fadeout  (S-SWF)  and  the 
sudden  cosmic-noise  absorption  (SCNA).  Furthermore,  the  increased 
ionisation. in. the  0-region  enhances  the  ionospheric  current  systems  and 
produces  a  discontinuity  in  the  magnetic  record  called  a  crochet.  These 
effects  last  only  as  long  as  the  flare  and  are  caused  by  an  increase  in  soft 
X-ray  radiation  (Chubb,  et  al.  ,1957;  Kreplih,  1961). 

Another  manifestation  of  solar  activity  is  in  the  radio  spectrum.  In  addi¬ 
tion  to  the  thermal  blackbody  radiation,  nonthermal  emissions  occur  at 
very  irregular  intervals  with  intensities  of  perhaps  10^  times  as  large  as 
‘Jtermal.  The  classification  of  radio  emissions  in  the. band  25  to  10,  000 
Vic  is  made  in  terms  of  the  following  spectral  types. 

Type  I  --  random  bursts  of  noise  which  may  last  for  hours  ''x  days  -- 
associated  with  sunspot  areas  but  not  particularly  with  flares. 

Type  II  —  slow-drift  burst  —  starts  at  high  frequencies,  drifts  stva.'Uy 
down  at  about  1  Me  with  a  duration  of  several  seconds. 

Type  III  —  fast  drift  dust  —  starts  at  high  frequency  with  a  duration  of 
several  seconds. 

Type  IV  —  continuum  radiation  --  strong  -  long  with  a  duration  of 
several  hours. 
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Type  V  —  continuum  .occurring  immediately  after  type  in  with  a  dura¬ 
tion  of  several  minutes. 

2.  Riometer  Data 

Solar  cosmic  rays  (relativistic  protons)  are  generated  in  significant  numbers 
at  the  time  of  large  solar  flares-.  They  reach  the  vicinity  of  the  Earth  within 
about  a  few  hours  after  the  beginning  of  the  flare.  The  trajectories  of  the  in¬ 
cident  protons  maybe  described  using  the  theory  of  Stormer  (1955)  ’..lich 
treats  them  as  individual  particles.  Because  of  the  guidance  of  the  magnetic 
field,  solar  cosmic  rays  are  confined  to  the  region  greater  than  60  degrees 
geomagnetic  latitude.  Production  of  ionized  particles  in  the  D-region  of  the 
ionosphere  causes  strong  absorption  of-HF  signals  reflected  from  the  iono¬ 
sphere.  Blackout  conditions  for  HF  radio  are  caused  in  the  polar  regions 
and  result  in  an  absence  of  ionospheric  reflections  in  high-latitude46nosonde 
recordings.  This  continues  for  about  1  day  and  is  termed  a  polar  cap  absorp¬ 
tion  event  (PCA  or  type  III  absorption). 

lonosonde  measurements  provide  only  a  qualitative  measurement  of  absorp¬ 
tion,  especially  as  the  blackout  condition  is  dependent  on  the  F2  critical  fre¬ 
quency,  and  in  any  event,  places  only  a  lower  bound  on  the  D-region  absorp¬ 
tion.  A  more  efficient  instrument  for  measuring  absorption'has  been  dubbed 
a  riometor  (1.  e. ,  relative  ionospheric  opacity  meter).  The  device  is  a  re¬ 
ceiver  which  continuously  recorde-hosmic  holsc'power.  In  the  Arctic,  fre¬ 
quencies  of  about  2? -Me  are  used  while  18  to  20  Me  are  usually  employed  at 
lower  latitudes..  At  these  frequencies  terrestrial  sources  are  quite  weak, 
and  the  cosmic  noise  sources  are  paramount.  This  extraterrestrial  noise 
can  penetrate  the  ionosphere,  since  the  measurement  is  made  above  the 
F-region  critical  frequency.  Under  normal  conditions,  the  nighttime  D-region 
absorption  is  negligible  and  provides  a  system  of  calibration.  Daytime 
photo-ionization  of  the  D  region  produces  absorption  and  reduces  the  galactic 
signal  power  which  penetrates  the  D  region.  The  noise  power  which  pene¬ 
trates  the  D  region  is  a  measure  of  the  integrated  absorption, 

Riometer  absorption  is  usually  quoted  in  decibels  per  pass  it  the  measuring 
frequency.  The  riometer  is  ure/ul  also  for  identifying  other  absorption 
events  in  the  auroral  zone  (type  II  absorption).  PCA  (typo  III)  is  smooth 
and  long  lasting  compared  to  type  II  which  is  associated  with  visual  aurora. 

3.  Magnetic  Data 

At  the  time  of  optical  flares,  giant  clouds  of  low-enerip^  gas  ar--  emitted 
from  the  sun.  Moving  through  the  region  between  the  sun  and  Earth  at 
about  1000  km/sec,  these  gas  clouds  reach  the  Earth  2  days  later,  on  the 
average  (Obayashi,  1961).  The  average  energy  of  the  protons,  neasured 
by  Explorer  X  at  about  22  earth  radii,  is  of  the  order  of  1  kev  (Rossi,  1962). 
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The  interaction  of  the  shock  wave  associated  with  the  gas  cloud,  the  cloud 
particles,  and  the  Earth's  nnagnetic  field  produces  a  geomagnetic  storm. 
Magnetic  records  serve  to  characterize  this  disturbance:  three  components 
are  measured,  but  the  horizontal  (H)  component  is  most  significantly  af¬ 
fected.  The  advancing  shock  front  causes  a  sudden  increase  in  the  H  com¬ 
ponent  -  the  sudden  commencement  (sc)  of  the  geomagnetic  storm.  For 
several  hours,  the  field  remains  above  normal’  (the  initial  phase)  plunges 
sharply  below  normal  (main  phase)  and  then  returns  slowly  to  its  prestorm 
condition  after  about  30  hours.  The  actual  magnitude  of  these  effects  is 
individual  with  each  storm  but  the  pattern  is  similar.  Tabulations  of  sud¬ 
den  commencements  and  other  magnetic  disturbances  are  made,  but  copies 
of  individual  station  magnetic  records  (magnetograms)  are  difficult  to  dis¬ 
tribute.  In  order  to  provide  workers  with  a  measure  of  the  disturbance 
variation,  several  indices  of  activity  have  been  devised. 

1 .  The  C-figure  is  an  arithmetic  mean  of  the  subjective  classification 
by  all  observatories  of  each  day's  magnetic  activity  on  a  scale  of  0 
(quiet)  to  2  (storm). 

2.  Kp  is  the  mean  K  index  from  12  observatories  between  47  to  63  de¬ 
grees  geomagnetic  latitude.  The  scale  ranges  from  0  (very  quiet)  to 

9  (extremely  disturbed).  It  is  a  measure  of  the  3-hour  range  of  the 
most  disturbed  component  of  the  field.  The  scale  is  adjusted  to  com¬ 
pensate  for  the  latitudiiuil  variation  of  the  disturbance. 

3.  Ap  is  a  daily  index  of  magnetic  activity.  It  is  computed  as  the  mean 
of  8  values  of  a  3-hour  index  Sp  which  is  1/2  the  average  gamma  range 
of  the  most  disturbed  of  the  three  magnetic  components.  In  practice 
this  is  computed  from  Kp  by  the  schedule: 

Kp  0123456  7  8  9 

A„  0  4  7  15  27  48  60  132  207  400 

P 

Kp  is  the  most  useful  for  ionospheric  studies. 
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IV.  HF  EQUIPMENT  AND  SYSTEMS 


A.  EQUIPMENT 

While  the  ionosphere  provides  a  reflecting  surface  enubltr.g  radio  signals  to  be 
transmitted  beyond  the  horizon  and  its  absence  precludes  all  but  scatter  prop¬ 
agation,  the  electronic  equipment  which  generates,  transmits,  receives,  and 
detects  is  also  essential.  Each  communication  system  possesses  certain 
characteristic  assets  and  liabilities  and  must  be  matched  to  the  ionospheric 
conditions,  transmission  path,  and  information  rate. 

Communication  ground-to-air  is  now  almost  exclusively  voice,  while  ground- 
to-ground  is  transmitted  in  CW,  AM,  DSB,  SSB,  and  FSK.  Several  communi¬ 
cations  systems  presently  in  use  are  tabulated  in  table  XI. 

B.  HF  SYSTEMS 

1,  U.S.  Coast  Guard 

In  the  USCG  HF  communication  system,  which  is  used  in  connection  with 
the  North  Atlantic  weather  patrol,  there  are  8  stations  -  ships  A  through 
E:  Argentia,  Newfoundland;  Bermuda;  and  Washington,  D.  C.  These  are 
shown  in  figure  39.  The  power  output  for  the  ships  is  0.  5  kw  and  for  shore, 
3  and  1 3  kw.  At  Washington,  a  monorhombic  antenna  is  used.  Normally, 
communication  from  weather- ships  is  on  approximately  4,  6,  8,  12,  and 
16  Me.  When  relays  go  through  Argentia  the  frequencies  used  are  approxi¬ 
mately  0,  5,  2,  4,  and  7  Me.  Two  frequencies  are  utilized  simultaneously 
to  avoid  cutoff.  During  disturbed  periods  it  is  the  practice  to  use  relays 
on  north-south  paths,  that  ic,  messages  from  station  A  go  through  C,  and 
messages  froni  B  go  through  Argentia,  Newfoundland,  or  through  station 
D  or  E.  Routing  through  Bermuda  is  done  only  when  propagation  condi¬ 
tions  require  it. 

2.  Norlant  -  1958 

The  HF  narrow  band  system  in  1953  serve-’  US  Air  Force  installations  in 
the  North  Atlantic  region  and  in  northern  parts  of  North  America.  Much 
of  the  communication  in  this  latter  region  was  in  connection  with  IGY  ac¬ 
tivities,  The  geographical  layout  of  this  system  is  shown  in  figure  40. 

The  numbers  beside  the  lines  indicate  the  number  of  assigned  frequencies 
for  that  link.  Many  of  these  stations  are  equipped  with  air-ground  and 
ground-air  equipment.  Information  on  this  type  of  equipment  may  be  found 
in  the  previous  chapter.  In  connection  with  networks,  only  point-to-point 
equipment  need  be  considered.  These  equipment  parameters  are  given 
in  table  VUI. 
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TABLE  VIII 


HF  NARROW-BAND  SYSTEM  INFORMATION 


Point-to-Point.  Sygtems 
Transmitter  LD-2 

Power  output;  up  to  90  kw  peak  envelope  SSB 
Frequencies:  10  preset  Irequency  select  any  of  six 
Functional  description:  transmit  four  3  kc  channels  simultaneously 
Antennas:  Z  Rtibmbics,  high  band  and  low  band,  the  installation  of 
which  is  designed  for  particular  paths, 

t 

Receiver 

(1)  R-390/URR 
Sensitivity:  0.  5jiv 

(2)  Hughes  EC-ISO 
Sensitivity;  0.  3jiv  <10  Me 

O.S^v  >10  Me 

Functional  description:  band  divided  into  four  3-kc  channels, 
Antennas;  Rhombics  in  pairs  for  high  and  low  bands  for  space 
diversity. 

(3)  SP-600 
Similar  to  R-390 
Approximate  sensitivity:  0,$nv 

Ground-Air  Systems 

Transmitter  MW-2  Westinghouse 

Power  output:  SSB  (upper  or  lower) 

2,  5  kw  peak  envelope 
AM  1.  25  kw  currier 
Frequency  range:  2, 0  to  30,  0  Me 

Antennas:  5  vertical  towers  ar  antennas  to  cover  5  frequency 

bands.  Towers  and  antennas  are  each  proportional  to 
wavelength. 

Backup:  LDT-2 
Power  output:  4  kw  peak 

Receiver 

Sensitivity;  0. 25  itv 

Functional  Description:  6  channel  --  each  channel  can  receive 

both  upper  and  lower  side  band  separately. 
Frequency  Range:  4.  0  to  30. 0  Me. 

Antennas:  4  Granger  log-periodic  curtai:  array,  vertical  polariza¬ 
tion,  1 10-degree  beamwidth  +  9  db  gain  over  isotropic. 
For  omnidirectional  coverage,  discone  antennas  are  used. 
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Figure  39  USCG  NORTH  ATLANTIC  WEATHER  PATROL.  HF 
COMMUNICATION  SYSTEM 
63-3576 
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Figure  40  NORLANT  •  1958  HF  COMMUNICATION  SYSTEMS 

63-3575 
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Figure  41  NORLANT  -  1961  HF  COMMUNICATION  SYSTEMS 
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In  practice,  this  netyork  is  used  together  with  others,  such  as  VHF  or 
land-line^networks.  .'^'or  this  one,  only  one  frequency  at  a  time  is  used:^so 
-that-,the-aosfr  of-operatj  inican-be-kept-to-aiminimum,  -For-the  -major  relay- 
station,  i.  e. ,  Harmon,  'Goose,  and  Lajes,  relaying  is  accomplished  elec¬ 
tronically.  If  one  link  is  saturated  then  another  is  selected, 

3.  Norlant  -  1961 

This  system  is  essentially  the  same  as  in  1958.  The  major  difference  is 
that  the  number  of  stations  has  been  somewhat  reduced,  This  network  is 
shown  in  figure  41.  The  equipment  installed  at  these  stations  is  about  the 
same  as  before.  Circuit-information,  such  as  the  approximate  assigned 
frequencies,  bandwidth  and  power  is  given  in  table  IX, 

Some  additional  comments  on  the  operation  of  this  network  are: 

a.  The  peak  load  on  common  user  traffic  is  from  about  1400  to  about 
2200,  local  time.  Channel  utilisation  is  about  50  percent  if  averaged 
over  the  whole  day.  If  a  channel  is  assigned  to  only  one  user  and  his 
utilization  exceeds  50  percent,  it  is  likely  that  during  peak  hours  a 
backlog  of  messages  will  develop. 

b.  The  average  length  of  a  message  is  about  ten  lines,  but  this  varies 
from  one  message  to  another.  Messages  involving  air  traffic  are 
usually  shorter  and  have  an  average  length  of  about  two  lines, 

c.  Relays  require  less  than  1  minute  or  so  handling  time.  Usually 
the  whole  message  is  received  before  retransmission.  Some  messages 
go  through  several  links  without  any  manual  handUng  -  reception  and 
retransmission  is  accomplished  automatically.- 

d.  Destination  of  traffic  is  highly  variable.  Many  agencies  use  this 
network  and  other  ones  which  are  connected.  So  traffic  flows  o"  a 
world-wide  scale,  with  some  messages  passing  through  many  links, 
say  6,  and  others  through  only  1, 

4,  Strategic  Air  Command  HF  Communications 

For  the  SAC  constant-alert  system  and  other  SAC  training  operations,  there 
are  4main  control  stations  in  the  United  States:  Westover,  Massachusetts, 
Barksdale,  Louisiana,  Omaha,  Nebraska,  and  March  Air  Force  Base, 
California.  These  are  interconnected  and  serve  to  control  SAC  aircraft. 
Aircraft  flying  in  northern  regions,  or  in  Europe  for  example,  receiv,  HF 
signals  directly  from  the  US  or  from  other  nearby  aircraft,  or,  if  possible, 
from  a  ground  station  on  UHF.  Sometimes,  subsidiary  stations,  such  as 
Goose  Bay,  transmit  HF  and  receive  and  relay  back  to  the  US.  The  US 
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TABLE  IX 


NORLANT  HF'CIRCUIT  INFORMATION  -  196 T 


System 

Allocated 

Frequencies 

(Me) 

Bandwidth 

(kc) 

Maximum 

Power 

(kw) 

Andrews  /Harmon 

5,  8,  14,  16. 

3 

3 

Harmon/Andrews 

5,  9.  12 

3 

3  to  4 

Andrews /Lajes 

6,  8,  10,  15,21 

3  to  12 

40 

Lajes /Andrews 

8,  10,  16,  19,22 

6  to  12 

4n 

Thule  /Sondre  Strom 

3,  5,  8,  9 

3 

3 

Sondrestrom/Thule 

5.  7,  8,  9 

3 

3 

Goose /Sondrestrom 

6,8,9,16 

3 

3 

Sondrestrom/Goose 

5,  7,  9,  12 

3 

3 

Harmon/Lajcs 

5, 8,  9, 10, 12, 16,  20 

3  to  9 

4  to  40 

Lajes /Harmon 

8,  10 

9 

10 

Lajes /Keflavik 

3,  4,  8,  12,  14,  20 

3 

3  to  10 

Kellavik/Lajes 

3,4,  9, 12,  16,  20 

3 

2  to  10 

Lajes /Croughton 

6,9,  16,  23 

3  to  12 

4  to  io 

Croughton/Lajes 

5,  8,9,  10,  14,  .16 

6  to  12 

4  to  40 

Lajes/Torrejon 

9,  20, 24 

6  to  12 

4  to  40 

Torrejon/Lajes 

5,  9,  10.  14,20,  24 

6  to  9 

3  to  40 

Lajes/ Sidi  Slimane 

5,  11, 12, 14, 16 

3 

3  to  5 

Slid!  Slimane/Lajes 

5,  7,9, 11,  14,  16 

3 

3 

i 
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stations  have  a  power  output  of  around  10  kw  and  up  to  45  kw,  if  necessary. 
Normally,  only  the  power  actually  required  is  used.  There  are  8  frequen¬ 
cies  available  in-therange-2-to-30  Me..  The-antennas-are  different-types: 

li  Disk  cones,  gain  2  db 

2,  Billboards,  gain  20  db 

3.  LPS  (rotating)  gain  13  db 

Aboard  the  aircraft  {B52,  KC135,  and  B47)  there  are  Collins  A  58  re¬ 
ceivers,  Their  transmitter  output  is  about  1-kw  peak  power.  They  use 
the  same  frequencies  as  the  control.  It  is  not  necessary  for  control  to 
monitor  more  than  four  frequencies  at  a  time,  because  of  the  diurnal  vari¬ 
ation  of  frequency  assignment.  In  actual  operation  it  would  be  preferable 
to  have  a  communication  operator  aboard,  but  the  position  has  been  aban¬ 
doned,  Thus,  the  present  crews  have  a  difficult  time  maintaining  communi¬ 
cation  while  flying  the  aircraft.  Besides  weather  information  and  mission 
data,  periodic  reports  on  SAC  aircraft  are  maintained  by  US  control.  In 
northern  regions,  use  is  made  of  DEW  line  facilities  for  UHF  contact. 

This  provides  a  geographical  band  for  high-flying  aircraft  about  250  miles 
on  either  side  of  the  line.  Further  north,  HF  is  the  only  rneans  of  com¬ 
munication.  Since  July  1961,  when  the  completely-operational  24-hour- 
air  alert  was  started,  there  have  been  no  serious  outages.  However,  there 
have  been  some,  and  they  result  from  propagation  conditions,  aircraft 
communications  equipment  failures,  and  nuclear  testing. 
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V.  AVAILABLE  CHANNELS  -  MAXIMUM  UTILIZATION 


A.  COMEU.TER-PROGRAM 

We  consider  a  communications  network  consisting  of  n  receisfing-transmitting 
stations,  each  represented  by  an  index  i,'i  ••  .  Station  i  is  capable  of 

transmitting  cjj  message-units  to  stations  j,  i,  j  •  1 . n,  during  the  period  be¬ 
ing  considered  and  it  is  desired  to  transmit  from  i  to  j,  i,  j  >  i  f  j, 

either  directly  or  through  intermediate  stations,  the  accumulated  number  of 
/,j  message-units. 

Two  related  problems  are  considered.  First,  a  determination  must  be  made  of 
the  largest  number  (or,  more  generally,  the  largets  value)  of  the  accumulated 
messages  that  can  be  transmitted  through  the  network  to  the  desired  destinations. 
Next,  having  determined  that  number  (or  value),  a  way  must  be  found  of  channel¬ 
ing  the  messages  through  intermediate  points,  when  necessary,  so  as  to  transmit 
the  information  at  a  minimum  cost. 

In  the  practical  cases  to  be  considered,  it  is  expected  that  many  of  the  values 
c;j  will  be  0,  indicating  that  station  i  has  no  capability  of  transmitting  to 
station  j  for  many  couples  (i,  j).  This  fact  permits  us  to  formulate  the  problem  in 
a  manner  that  will  tend  to  render  our  computations  more  manageable. 

Assuming  that  an  ordered  couple  (i,  j)  of  stations  is^connected  by  a  link,  if 
Cjj  -  0 --There  may  be  two  district  links,  (i,  j)  aiwi  0i  •)  with  cjj  4  cjj  ,  in  general; 
or  (i,  j)  may  be  connected  by  a  link,  but  not  (j,  i)--  the  links  may  be  enumerated  in 
some  arbitrary  fashion  and'represent  them  by  an  index  I,  I  -  I,..., I,  letting  Cf  »  cjj, 
for  the  link  /  which  connects  i  to).  Furthermore,  the  addresses  of  messrges 
shall  be  enumerated,  each  address  representing  an  ordered  couple  (i,  j)  such 
that  /jj  i  0  thus,  an  address  represents  an  origin-destination  combination 
characterizing  a  given  message  or  a  combination  of  a  forwarding  and  a  return 

address.  Let  /„  a  »  1 . a,  where  /,  is  the  nmnber  of  message  un’ts  to  be 

sent  according  to  the  address  a. 

Finally,  the  various  admissible  paths  which  a  message  may  take  wlien  bearing 
an  address  a  must  be  described.  A  gi/en  path  mr.y  he  uniquely  defined  by 
indicating  the  links  of  which  it  is  composed.  The  paths  are  enumerated  by  an 
index  p,  p  °  1,— ,  p,  grouping  the  paths  belonging  to  a  given  address  next  to 
each  other  and  enumerating  those  pertaining  to  address  a  +  1  right  sifter  those 
pertaining  to  address  a,  a  »  1,  2 . a  -  1. 

A  message  transmitted  through  a  path  p  may  have  a  given  value.  In  general, 
one  would  tend  to  make  all  values  equal  unless  certain  addresses  are  to  be 
favored  (i.  e. ,  for  priority  reasons)  or  even  certain  paths  are  to  be  favored 

(say,  for  strategic  reasons).  Let  v^,  p  °  1 . p  be  the  value  of  path  p,  which, 

unless  otherwise  indicated,  is  set  equal  to  1. 
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Different  costs  of  transmitting  a  message  through  a  given  path  p  may  be  con¬ 
sidered.  Thus,  as  an  example,  the  cost  of  transmitting  a  message  from  one 
station  to  another  may  vary,  yielding  different  costs  for  different  links;  in 
that  case,  the  cost  of  transmitting  a  message-unit  through  path  p  is  the  sum  of 
the  unit-costs  associated  v/ith  the  links  of  which  the  path  is  made  up.  Unless 
otherwise  specified,  let  the  unit-cost  for  path  p  be  equal  to  the  number  of  links 
in  it,  making  the  cost  proportional  to  the  labor-hours  required  to  transmit  a 
message-unit  through  the  path.  The  unit  cost  for  path  is  represented  by  Wp. 

1  "  1 . P 


Now,  Xp.  p  <•  1.—  p>  represent  the  number  of  message-units  sent  through  the 
path  .  Let  a/  p  -  1  if  link  I  belong’  to  path  P  and  let  otherwise  p  -  0  .  Let 
a^  1  1  if  the  path  p  belongs  to  address  a  and  let  otherwise  p-O* 

The  first  problem  is  that  of  determining  the  maximum  value  of  accumulated 
messages  that  can  be  transmitted  through  the  network:  the  problem  can  be 
formulated  asrfollows:  find  nonnegative  Xp,  p  •  i....,p,  such  that  the  function 

p 

V  >  Sp.  j  Tp  Xp  is  maximxim  subject  to  the  restrictions; 


•/,  p  *p 


<  C^,  /  ■  It  ...t  / 


(64) 


and 


\  «  1  *7  +  a,  p  *p  -  ^a'  *  ■  •••• 


(65) 


The  quantity  v  represents  the  total  value  (or,  in  a  special  case,  the  total  number) 
of  the  messages  transmitted  through  the  network.  Restrictions  (64)  specify  that 
the  tots.1  number- of  message  units  channeled  through  link  I  (which  may  be  shared 
by  many  paths)  cannot  exceed  cf .  Similarly,  restrictions  (65)  specify  that  no 
more  than/,  messages -units  are  available  for  useful  transmission  throv^h  the 
various  paths  servicing  the  address  «. 


After  this  first  problem  is  solved,  there  is  the  quantity  available  which  is 
the  maximum  value  of  v  subject  to  the  restticiions.  Nov/  solve  the  optimal 

channeling  problem  by  minimizing  the  total  cost  w  =  S  ^  ^  Wp  jl  subject  to  Xp 
being  nonnegative,  subject  to  restrictirns  (64)  and  (65)^and  subject  to  the  added 
restriction 


P-1 


X_ 

p  p 


''max 


(66) 


This  last  restriction  specifies  that  costs  may  not  be  cut  by  transmitting  less 
than  a  total  value  of  messages. 
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Both  of  these  problems,  as  formulated,  belong  to  the  class  of  linear  programming 
problems.  Illustrative  examples  were  solved  at  Avco  RAD  using  a  digital  com* 
puter  program  (Math  Section  program  1332)  based  on  a  convex  programming 
method  described  in  RAD-TM-62-58, 

B.  EXAMPLES 

To  show  how  the  mathematical  program  of  optimized  routing  is  carried  out, 
some  definitions  will  be  reviewed  and  then  three  examples  will  be  given.  The 
network  used  is  one  with  nine  nodes  or  stations  which  are  interconnected  as 
shown  in  figure  42.  Definitions  are  given  below  for  terms  which  appear  often 
in  the  test; 

1.  A  link  is  a  connection  between  two  stations.  The  links  for  the  example 
network  are  indicated  as  shown  by  the  lines  in  figure  42. 

2.  An  address  is  a  designation  of  the  source  and  destination  of  a  message. 

3.  A  path  is  the  route  by  which  a  message  is  sent.  In  general,  a  message 
will  be  sent  over  one  or  more  links. 

4.  The  load,  Ijj.  is  made  up  of  all  the  messages  to  be  sent  to  an  address 
in  a  given  interval  of  time. 

5.  The  capacity,  cjj,  is  made  up  of  all  the  messages  which  can  be  sent 
over  a  link. 

In  HF  communication  problems  it  will  be  ij-  and  cj^  which  are  used  as  inputs  to 
an  optimization  program.  Ijj  will  be  determined  by  the  users  of  the  communica¬ 
tion  facility,  and  C|^  will  be  determined  by  the  equipment  being  used  and  by  the 
condition  of  the  ionosphere.  In  these  examples,  both  the  i..',  and  the  will’ 
be  assumed  known.  For  the  first  example  the  value  of  C|^  is  taken  to  be  10  for 
all  24  links  and  the  1.^'$  are  listed  in  table  X.  The  total  is  34  messages  which 
may  be  interpreted  as  a  message  rate,  or  .as  pages  per  hour. 
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TABLE  X 


Address 

Messages 

hi 

‘1} 

1.7 

4 

3.5 

7 

4.1 

6 

.  4,5 

3 

6.7 

5 

6.8 

9 

Total 

34 

For  these  examples  the  units  are  not  so  important  as  the  relative  values  of  ijj 
and  This  example  is  a  case  when  the  capacity  of  the  network  is  somewhat 

larger  than  the  load  placed  on  it.  It  would  be  expected,  then,  that  the  paths 
chosen  would  be  the  most  direct  and  that  they  could  possibly  have  been  selected, 
by  inspection.  The  routing  of  messages  from  6  to  7  and  from  6  to  8  is  shoy^n  in 
figure  43a,  as  a  82unple  of  the  results. 

The  routing  of  6,  7  is  obvious.  For  6,  8  it  is  evident  that  only  five  message 
units  can  pass  over  link  6,  7  and  that  the  remaining  four  messages  must  pass 
over  other  links.  Link  6,  7  is  the  only  one  which  becomes  saturated.  Since 
the  capacity  is  large  the  selection  is  somewhat  arbitrary.  All  four  messages 
could  just  as  well  have  been  sent  over  either  route  actually  selected  by  the 
computer.  If  these  links  are  open  for  traffic  to  pass,  it  should  make  little 
difference  what  is  the  choice.  However,  if  it  were  desirable  to  eliminate 
this  arbitrariness,  criteria  could  easily  be  set  up  to  pick  a  preferred  route. 

The  other  messages  to  be  sent'in  this  example  are  straightforward  and  there 
is  no  saturation  of  any  link  other  than  6,  7. 

In  the  second  example  the  link  capacity  has  been  kept  at  10,  but  the  number  of 
addresses  is  increased  to  12.  This  corresponds  to  a  case  when  the  capacity  of 
the  network  is  not  large  compared  to  the  load.  The  first  six  addresses  arc  -.s 
before.  These  Ijj  and  the  new  ones  are  shown  in  table  XI. 
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Figure  42  LINKS  FOR  THE  EXAMPLE  NETWORK 


Figure  43  ROUTING  OF  MESSAGES  -  SAMPLE  RESULTS 


TABLE  XI 


Address 

Messages 

Address 

Messages 

i.  i 

1 

i 

1 

1.7 

4 

7,1 

2 

3,5 

7 

7,2 

4 

4.1 

6 

8,1 

6 

4.5 

3 

8,5 

1 

6.7 

5 

8.9 

8 

6,8 

9 

9,4 

2 

ToUl  1, 

j  57 

One  important  difference  between  the  results  of  this  example  and  those  of  the 
first  is  that  four  message  units  could  not  be  sent  out  of  the  total  of  57.  If  one 
looks  at  the  message  units  going  to  station  1  the  reason  becomes  obvious.  Out 
of  a  total  of  14,  only  10  can  pass  over  link  2,  1.  Tims,  four  of  the  message 
units  could  have  reached  only  as  far  as  station  2. 

In  other  respects  this  example  is  quite  similar  to  the  first  one.  The  results  for 
the  two  addresses  discussed  before  (6,  7  and  6, 8)  are  shown  in  figure  43b. 

The  routing  of  6,  7  is  again  obvious.  The  routing  of  6, 8  is  different,  but  this  is 
still  of  cm  arbitrary  nature.  The  capacity  after  other  messages  have  been  sent 
is  large  enough  on  path  6,  9,  7,  8  so  that  the  routing  could  have  been  the  same 
as  before. 

In  the  third  example,'  the  number  of  addresses  and  their  loads  has  been  kept  the 
same  as  in  the  second  example  (see  table  XIV).  However,  the  link  capacity  is 
reduced  to  four.  This  corresponds  to  a  case  when  the  capacity  of  the  neh/ork  is 
somewhat  less  than  the  load.  Now  the  routing  is  rather  sophisticated  and  could 
not  easily  be  anticipated  without  a  mathematical  program  such  as  the  one  used. 
Although  the  network  is  used  in  its  most  efficient  way,  23  message  units  could 
not  be  sent.  The  reason  is  similar  to  the  previous  case ,  but  is  not  easily  seen 
in  every  instance.  The  routing  for  the  two  addresses,  6,  7  and  6,  8  is  shown 
again  (figure  43c)  for  comparison  with  the  previous  examples. 

Here  the  routing  is  not  arbitrary  to  any  extent.  For  the  address  6,  7  link  6,  7 
becomes  utilized  fully,  so  one  message  unit  is  rerouted  --  the  shortest  way. 

For  the  address  6,  8  only  three  of  the  mne  message  units  could  be  sent.  The 
ones  sent  are  divided  according  to  the  available  capacity  on  the  paths  chosen. 

The  significemt  fact  relating  to  the  choice  of  paths  is  that  whenever  the  number 
of  messages  to  be  sent  can  be  increased  by  rerouting,  the  computer  will  indicate 
which  ones  are  to  be  used.  Itis  possible  to  spell  out  just  what  choices  are 
arbitrary  and  which  ones  are  not. 
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VI.  COMPLETED  SYSTEM 


To  show  how  the  concepts  and  Infornnation  described  in  this  report  utillsed»  an 
example  of  a  hypotheticaL  but  realistic  commxmi cation  ni'oblerti  will  be  given. 
Firsti  the  network  and  equipment  parameters  will  be  assumed,  Then>  the 
available  frequencies  will  be  found  as  would  exist  during  a  typical  polar-cap 
absorption  event.  Finally,  for  a  given  message  load,  the  optimum  routing  and 
selection  of  frequencies  will  be  carried  out. 

The  network  used  in  this  example  is  the  same  one  used  in  the  previous  chapter, 
but  it  has  been  superimposed  on  a  map,  (figure  44)  so  that  realistic  communica¬ 
tion  facilities  can  be  imagined.  Assignment  of  equipment  to  these  stations  is 
expressed  in  terms  of  equipment  parameters  rather  than  their  names  or  types. 
These  parameters  are  listed  in  table  XU.  From  these  values  it  is  possible  to 
compute  the  radiated  field  strength  at  1  km  from  the  transmitter,  Eg  ,  in  the 
direction  of  the  receiver.  The  transmitter  antennas,  receiver  antennas,  and 
receiver  sensitivity  at  a  station  are  assumed  the  same  for  all  links  in  and  out 
of  that  station,  so  that  this  information  can  be  listed  by  station  rather  than  by 
link.  In  fact,  tltese  parameters  have  been  chosen  to  be  the  same  at  all  stations, 
but  there  is  no  particular  advantage  in  doing  this.  Eg  .  listed  in  table  XIII,  is 
given  by  74. 8  +  10  log  Pg  where  Pg  »  Pt/3  nnd  P,  is  the  power  output  in  watts. 
The  other  parameter  listed  in  table  XIU  is  the  minimum'  signal  level  receivable 
if  the  signal-to-noise  ratio  is  to  be  10.  These  also  can  be  listed  by  station, 
but  normally  they  would  be  listed  by  link.  The  allowable  loss  of  signal  strength, 
for  any  mode  which  is  possible  by  consideratiou  of  the  reflection  layers,  is 
then  determined  by  the  difference  between  the  signal  level  indicated  by  Eg  at 
the  transmitting  site  and  the  signal  level  Indicated  by  the  minimum  sigMl  level 
at  the  receiver.  This  difference  must  be  adjusted  to  take  into  account  ^e  re¬ 
ceiver  antenna.  The  allowable  loss  is  found  in  column  three  of  table  XIV.  The 
adjustment  F  is  found  from  F>  30.  6  -  20  log  f  ,  with  {  in  Me.  The  space  at¬ 
tenuation  is  also  given  here.  When  that  ie  deducted  from  the  allowable  total 
loss,  the  allowable  loss  due  to  absorption,  reflection,  and  ground  losses  is 
found.  In  this  example,  the  latter  two  are  neglected  for  simplicity.  It  might 
be  appropriately  mentioned  here  that  sporadic  E  modes  are  also  nei,lected. 

As  shown  earlier,  it  is  in  relation  to  the  presence  of  sporadic  E  that  reflection 
losses  become  significant.  So,  for  this  example,  the  allowable  loss  due  to 
absorption  alone  is  found  after  space  attenuation  is  deducted.  This  absorption 
loss  is  also  included  in  table  XI'V. 

A  number  of  frequencies  are  assigned  to  each  link.  The  assignment  of  these 
frequencies  was  guided  by  some  general  consideration,  but  the  details  were 
more  or  less  arbitrary.  The  lowest  frequency  assigned  is  3  Me  with  the  highest 
around  25  Me.  The  spacing  between  frequencies  for  any  one  link  was  about  4 
or  5  Me.  For  each  reciprocal  link  the  assigned  frequencies  differed  by  0. 1  Me 
in  every  case.  It  will  also  be  noticed  that  for  the  shortest  links,  the  low 


-107- 


TABLE  XII 


SYSTEM  PARAMETERS  OF  A  HYPOTHETICAL 
COMMUNICATION  NETI^ORK, 


Station 

Pt 

il 

R.S. 

Gr 

Total  Noise 

watts 

El 

/iV 

db 

db  above  1  nr 

3x  10^ 

10 

0,5 

0 

0 

3 

10 

0.5 

0 

0 

3 

10 

0.5 

0 

10 

3 

10 

0.5 

0 

10 

10 

10 

0..5 

0 

15 

3 

10 

0.5 

0 

5 

0.  5 

10 

0.5 

0 

-10 

8 

0.5 

10 

0.5 

0 

-10 

9 

3. 

10 

0.5 

0 

10 

P,  Tran*mltter  Power 

G,  Transmitter  Antenna  Gain 

G]^  Receiver  Antenna  Gain  over  a  half-wave  dipole 

R.S.  Receiver  Sensitivity 


TABLE  Xm 

RADIATED  FIELD  STRENGTH  AND 
MINIMUM  SIGNAL  LEVELS 


Station 

pQ 

db  above  1  fir  fm 

Minimum 
Signal  Level 
db  above  1  in 

1 

V.t 

10 

2 

116 

10 

3 

116 

20 

4 

116 

20 

5 

121 

25 

6 

116 

15 

7 

4 

8 

108 

4 

9 

116 

20 
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TABLE  XIV 

ALLOWABLE  ABSORPTION 


Eo 

Minimum 

Signal 

Level 

Allowable 

Lose 

Space 

Attenuation 

Allowable 

Absorption 

mm 

M 

51+F 

bSI 

IRI 

51 

HSH 

116 

59 

37 

2  to  5 

116 

65 

26 

2  to  8 

116 

112 

60 

52 

2  to  9 

116 

20 

96 

63 

33 

3  to  2 

116 

10 

106 

59 

47 

3  to  4 

116 

20 

96 

60 

36 

4  to  3 

116 

20 

96 

60 

36 

4  to  5 

116 

25 

91 

59 

32 

5  to  2 

121 

10 

111 

65 

46 

5  to  4 

121 

20 

101 

59 

42 

5  to  6 

121 

15 

106 

59 

47 

6  to  5 

116 

25 

91 

59 

32 

6  to  7 

116 

112 

54 

58 

6  to  9 

116 

96 

62 

34 

7  to  6 

108 

93 

54 

39 

7  to  8 

108 

104 

55 

49 

7  to  9 

108 

88 

59 

29 

8  to  2 

108 

98 

60 

38 

8  to  7 

108 

104 

55 

49 

9  to  2 

116 

106 

63 

43 

9  to  6 

116 

101 

62 

39 

9  to  7 

1 _ 

116 

mM 

112 

59 

1 

” _ 1 

frequencies  were  stressed  and  for  the  long  links  the  high  frequencies.  However, 
for  the  long  links,  low  frequencies  were  also  assigned.  A  complete  list  of  the 
assigned  frequencies  is  provided  in  table  XV.  It  may  be  noticed  that  there  is 
some  similarity  between  these  frequencies  and  the  ones  assigned  in  the  Norlant 
circuit  information  given  in  section  IV. 

The  next  step  is  to  find  which  of  these  frequencies  are  available.  A  single  set 
of  ionospheric  maps  has  been  chosen  for  this  example.  These  maps  are  for 
1500  UT,  11  April  1958  when  a  moderate  to  strong  PCA  event  was  taking  place. 

It  is  assumed  that  the  conditions  represented  on  the  maps  (figure  45)  are  valid 
for  a  certain  time  interval,  say  1  hour  (a  half  hour  before  and  c  half  hour  after 
the  map  time).  The  available  frequencies  found  during  this  ionospheric  situa¬ 
tion  will  determine  the  network  capacity  for  this  time  interval.  These. fre¬ 
quencies  are  found  in  two  steps. 

In  the  first  step,  analysis  is  made  of  the  range  of  frequencies  which  are  possible 
according  to  the  reflecting  properties  of  the  ionosphere  without  regard  to  prop¬ 
agation  losses.  In  this  study,  these  ranges  are  solely  determined  by  ray-path 
geometry  and  the  critical  frequencies  of  the  various  layers.  The  methods 
employed  are  described  in  detail  in  section  II.  Calculation  of  these  frequency 
ranges  yields  the  results  shown  in  table  XVI. 

For  the  simple  network  used  in  this  example  only  a  few  modes  are  available. 
These  are  the  single  E,  Fl,  and  F2  hops,  and  the  double  and  triple  Ehop.  As 
expected  the  lowest  frequencies  are  possible  for  E  hops  (Modes  1,  4,  and  8), 
a  narrow  frequency  range  for  Fl  hops  (Mode  2),  and  the  highest  frequenj^es 
for  F  hops  (Mode  3), 

In  the  second  step,  propagation  losses  are  calculated  for  the  assigned  frequencies 
that  lie  within  the  frequency  range  available  for  each  mode.  Calculation  is 
carried  out  using  the  methods  described  in  section  II.  The  results  of  those 
calculations  are  shown  in  table  XVll.  The  absorption  is  listed  for  all  assigned 
frequencies  which  are  possible  according  to  I'ay-path  geometry.  V'hen  the 
absorption  is  greater  than  the  amount  allowable,  then  that  frequency  is  con¬ 
sidered  not  available.  No  assigned  frequencies  fall  within  the  range  of  possible 
frequencies  for  Fl  hops,  so  mode  2  does  not  appear  in  table  XVII.  The  values 
listed  under  Received  Field  Strength  Ui's  in  units  of  db  above  the  minimum  sig¬ 
nal  level. 

If  this  value  of  Received  Field  Strength  is  zero  then  the  signal-to-noise  ratio  is 
10,  so  there  is  no  excess  in  field  strength  over  the  minimum  required  tovel. 
Whenever  the  calculated  Received  Field  Strength  was  found  to  be  below  10,  the 
corresponding  frequency  was  considered  not  to  be  available  although  in  theory 
it  would  be.  This  procedure  allows  a  margin  of  error  in  computed  versus 
actual  values.  So  the  frequencies  available  are  the  ones  listed  in  table  XVni 
which  have  a  Received  Field  Strength  greater  than  or  eqvial  to  10  db  aoove  the 
minimum  signal  level.  These  frequencies  are  listed  in  table  XVIII  without 
regard  for  the  mode  applicable. 
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Figure  45  IONOSPHERIC  MAPS  FOR  1500  UT,  n  APRIL  1958 
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Figure  45  (Concl'd) 


FREQUENCY  RANGE  AVAll^BEE  WIIHOVT  PROPAGATION  LOSSES 


TABLtXVU 


COMI>UTED  VALUES  Of  AECRVED  FIELD  STAENCTH 


115 


>;SAT  'FREQUENCIES 

(Received  Signal  Strength  >  10  db  Above  Minimum) 

Link 

Usable  Frequencies 

Link 

Usable  Frequencies 

I  to  2 

7.0  12.0 

5  to  6 

15.0  19.2  23.2 

2  to  1 

7. 1  12. 1 

6  to  5 

15.  1  19.  3  23.  3 

2  to  3 

14.0  17.0 

6  to  7 

7.4  10.2  14.2 

2  to  5 

6  to  9 

20.4 

2  to  8 

18.0 

7  to  6 

14.  3 

2  to  9 

19.0 

7  to  8 

12.6  16.4 

3  to  2 

9.1  14.1  17.1 

7  to  9 

19.4 

3  to  4 

16.0  21.0  25.0 

8  to  2 

18.1 

4  to  3 

16.1  21.1  25.1 

8  to  7 

12.7  16.5 

4  to  , 5 

18.2  23.0 

9  to  2 

19.1 

5  to  2 

20.  5 

5  tc  4 

18.3  23.1 

15.3  19.5 

To  find  the  amount  of  material  which  can  be  sent  over  each  link  certain  assump¬ 
tions  are  made  regarding  the  communication  facilities.  For  each  link  a  number 
of  channels  are  said  to  exist.  If  each  of  these  has  a  capacity  of  60  words  per 
minute,  then  the  total  capacity  on  a  link  for  1  hour  is  the  number  of  channels 
times  3600.  In  table  XIX  the  channels  are  listed  in  column  one  and  the  link 
capacity  for  these  channels  expressed  in  thousand-word  units  are  listed  in 
column  two.  In  the  next  two  columns,  the  channels  and  capacity  available  are 
listed. 

The  assigned  load  to  be  sent  over  the  network  is  listed  in  table  XX.  Comp..'a- 
tion  can  then  be  made  on  how  the  load  should  be  distributed  for  the  available 
capacity.  This  can  also  be  done  for  the  case  when  all  the  channels  are  open. 

The  overall  results  can  be  expressed  in  terms  of  the  load  sent  or  the  load  re¬ 
maining!  the  latter  is  shown  in  table  XX  for  both  maximum  and  available  capac¬ 
ity. 

It  is  seen  that  the  message  load  which  can  be  sent  during  the  PCA  is  less  than 
the  amount  which  could  be  sent  if  the  full  capacity  of  the  network  were  available. 
This  condition  is  not  surprising,  as  it  occurs  at  certain  times  of  the  day  and 
year  even  in  undisturbed  conditions.  To  utilise  the  available  capacity  to  the 
best  advantage,  message  sent  of  the  highest  priority  should  be  sent  first. 

In  the  design  of  this  network  there  were  only  a  few  links  corrected  from  any  one 
station.  In  more  flexible  networks  other  links  which  are  found  to  offer  high 
capacity  can  be  set  up.  For  example,  equipment  could  be  installed  so  that 
direct  communication  can  be  made  between  stations  nine  and  four,  or  between 
six  and  two.  On  these  links,  there  may  be  other  modes  than  the  slniple  ones, 
such  as  M  hops  which  occur  often  during  the  critical  stages  of  ionospheric 
storms.  Only  by  extensive  testing  of  link  capacity,  or  frequencies  available, 
under  a  variety  of  storm  conditions  can  the  best  links  be  found.  Of  the  best 
ones  found,  those  most  consistent  with  actual  requirements  are  chosen.  Thus, 
as  would  be  expected,  the  final  design  of  a  network  depends  upon  the  expected 
loads  and  the  capacities  found  for  quiet  and  disturbed  ionospheric  conditions. 


TABLE  XIX 


LINK  CAPACITIES 


Link 

Maximum  Number 
of  16- channel 
Teletype  Connections 

Maximum 

Capacity 

•  1 

Available  Number 
of  16-channel 
Teletype  Connections* 

Available 

Capacity 

2 

115 

2 

115 

nin 

2  I 

115 

2 

115 

3  1 

173 

115 

2  to  5 

2 

115 

0 

0 

2  to  8 

1 

58 

1 

58 

2  to  9 

2 

115 

1 

58 

3  to  2 

2 

173 

2 

115 

3  to  4 

3 

173 

3 

173 

4  to  3 

3 

173 

3 

173 

4  to  5 

3 

173 

2 

115 

5  to  2 

2 

115 

0 

0 

5  to  4 

3 

173 

2 

115 

5  to  6 

3 

173 

3 

173 

6  to  5 

2 

173 

2 

115 

6  to  7 

2 

115 

2 

115 

6  to  9 

3 

173 

1 

58 

7  to  6 

2 

115 

1 

58 

7  to  8 

1 

58 

1 

58 

7  to  9 

2 

115 

1 

58 

8  to  2 

1 

58 

1 

58 

8  to  7 

1 

58 

1 

58 

9  to  2 

2 

115 

1 

58 

9  to  6 

3 

173 

1 

58 

9  to  7 

2 

1 _ 

115 

2 

115 

’I'This  number  is  equal  to  the  number  of  usable  frequencies  (see  table  XVIIl) 
or  to  the  maximum  number  of  I6>channel  teletype  connections,  whiche/er 
is  lowest. 
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